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• The hard-scattering contributions to heavy-to-hght form factors at large recoil are stud- 
ied systematically in soft-collinear effective theory (SCET). Large logarithms arising 

\ from multiple energy scales are resummed by matching QCD onto SCET in two stages 

^ ■ via an intermediate effective theory. Anomalous dimensions in the intermediate the- 



ory are computed, and their form is shown to be constrained by conformal symmetry. 
Renormalization-group evolution equations are solved to give a complete leading-order 
analysis of the hard-scattering contributions, in which all single and double logarithms 
are resummed. In two cases, spin-symmetry relations for the soft-overlap contributions 
to form factors are shown not to be broken at any order in perturbation theory by hard- 
scattering corrections. One-loop matching calculations in the two effective theories are 
performed in sample cases, for which the relative importance of renormalization-group 
evolution and matching corrections is investigated. The asymptotic behavior of Sudakov 
logarithms appearing in the coefficient functions of the soft-overlap and hard-scattering 
contributions to form factors is analyzed. 



1 Introduction 



Weak-interaction form factors for exclusive heavy-to-light transitions at large recoil energy, 
such as B ^ TT I ly with E^^ ~ itib/^, are an important input to measurements of the parameters 
of the unitarity triangle. The QCD description in this energy regime is complicated by the 
competition between different scattering mechanisms and the resulting proliferation of relevant 
energy scales. The tools of effective field theory provide an efficient means of separating the 
contributions from different scales and setting up a controlled expansion in the small ratios 
Aqcb/E and AqcoMfe- 

The appropriate theory in the present case is soft-coUinear effective theory (SCET), which 
is constructed to describe processes with both soft and coUinear partons [1, 2, 3, 4, 5]. Using 
SCET, it has been argued that there are two competing contributions to large-recoil heavy- 
to- light form factors at leading power in Aqqd/E (ignoring scaling violations from Sudakov 
logarithms) , referred to as the soft-overlap (or Feynman) mechanism and the hard-scattering 
(or hard spectator-scattering) mechanism [6, 7, 8]. In the first of these, the spectator quark 
in the S-meson is absorbed into the light final-state meson with no large momentum transfer. 
For this to happen, both the initial- and final-state partons must be arranged in an endpoint 
configuration with atypically small values of certain momentum components. In the second 
mechanism, a large momentum is transferred to the spectator quark via hard gluon exchange. 
The suppression due to the wave-function fall-off in the first case, and the suppression due to 
hard momentum transfer in the second case, are of the same order in power counting. 

In this paper, we present a renormalization-group (RG) analysis of the hard-scattering 
contributions. In SCET this mechanism is described by non-local four-quark operators, whose 
matrix elements factorize into products of leading-twist light-cone distribution amplitudes 
(LCDAs) for the S-meson and the light final-state meson. The matrix elements are multiplied 
by calculable coefficient functions, and the resulting convolution integrals are convergent to 
all orders in perturbation theory [7, 8] . The coefficient functions at an appropriate low-energy 
hadronic scale may be computed to any order in ag by perturbative matching of QCD onto 
the effective theory and subsequent RG evolution down to the low-energy scale. The analysis 
applies also to more comphcated decay processes such as S — > tttt and B — > K*^, for which 
QCD factorization formulae relate the decay amplitudes to the S ^ tt or S — > K* form 
factors plus a residual hard-scattering term [9, 10, 11]. 

The soft-overlap contributions to heavy-to-light form factors can be defined in SCET in 
terms of matrix elements of effective-theory operators obeying spin-symmetry relations ap- 
propriate for a heavy-coUinear transition current [8]. The relevant operators are rather com- 
plicated and lead to "non-factorizable" matrix elements sensitive to endpoint momentum 
configurations, transverse momentum components, and non-valence Fock states. However, 
these soft overlap contributions can be described in terms of universal functions C,m{E) that 
only depend on the light final-state meson but not on the Lorcntz structure of the currents 
whose matrix elements define the various form factors. For instance, there is one function 
(p{E) for decays into pseudoscalar mesons, and similarly only one function each, Cv||(-E') and 
(vj_{E), for decays into longitudinally and transversely polarized vector mesons. This implies 
spin-symmetry relations between different soft form factors, which were first derived in [12] 
by considering the large-energy limit of QCD. 
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The presence of both a soft-overlap and a hard-scattering contribution is summarized by 
the factorization formula [13] 

Ff-^^{E) = Q{E)CM{E)+ / — 0bM / dufM<pM{u)Ti{E,u;,u), (1) 

Jo ^ Jo 

which is vahd at leading power in Aqqd/E. Here 0b and 0m are the LCDAs of the B- 
meson and light final-state meson, and /m is the decay constant of the light meson. The 
Wilson coefficients of the effective-theory operators, Ci{E) and Ti{E,uj,u), may be calculated 
in perturbation theory to any order in as, and a RG analysis can be used to relate the 
coefficients at different renormahzation scales. However, the large- behavior of the soft- 
overlap contribution cannot be addressed satisfactorily in perturbation theory, since the long- 
distance matrix elements, (m{E), depend on the energy £^ in a non-pcrturbativc way [8]. Thus 
the issue of whether one of the soft-overlap or the hard-scattering contributions is enhanced 
relative to the other in the formal limit E ^ oo cannot be addressed using short- distance 
methods. Phenomenologically, it appears that the soft-overlap terms are dominant for physical 
values of the coupling and mass parameters. Although not a complete answer to the question, 
the relative suppression of the coefficients multiplying the long-distance matrix elements may 
be computed using perturbative methods. Studying the resummation of Sudakov logarithms 
for these coefficients, we find that the soft-overlap contribution is suppressed in the formal 
asymptotic limit E ^ oo, but that this suppression is mild for realistic values of E. 

Because the form factors involve three different physical scales, namely //^ ~ ml (hard), 
/i^ ~ mbAqcD (hard-coUinear) , and /i^ ~ ^qcd (soft), integrating out modes of progressively 
smaller virtuality results in a sequence of effective theories. At the high scale /x^ ~ mf, the 
effective theory is described by the usual QCD Lagrangian (with five quark flavors) plus an 
effective weak-interaction Lagrangian obtained by integrating out virtual W and Z bosons and 
top quarks. Integrating out modes of virtuality ml we arrive at an intermediate effective theory 
containing soft modes and hard-coUinear modes of virtuality pf^^. ~ mftAqcD- In this paper we 
are mainly concerned with this intermediate theory, called SCETj [2, 3, 4]. Integrating out the 
hard-collinear modes of virtuality mfoAqcD yields the final low-energy theory, denoted SCETn, 
consisting of soft and coUinear modes of virtuality p^, pi ~ Aq^j^. In this case, soft-coUinear 
messenger modes are also required [14, 15]. 

In the following section we briefly review some relevant elements of SCET. Section 3 lists 
the leading and subleading SCETj current operators, which are required for the discussion 
of heavy-to-light form factors, together with the matching coefficients for these operators 
obtained at tree level. For the example of the scalar current, we also calculate the one- 
loop matching coefficients of the leading and subleading operators. Section 4 contains the 
main result of the paper, i.e., the anomalous dimensions of the subleading SCETi current 
operators. These quantities exhibit an interesting symmetry property, which can be traced 
back to a conformal symmetry (at the classical level) in the hard-collinear sector. We briefly 
review the constraints imposed by conformal symmetry in Section 5, and use these results to 
diagonalize the non-local part of the evolution operator. We also present a formal algebraic 
solution of the evolution equations for the currents and their Wilson coefficient functions. In 
Section 6, we discuss the operator representation and renormalization for the hard-scattering 
contributions in SCETn. In two cases, namely for the form-factor ratios Ai/V and T2/T1, 
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we find that the spin-symmetry relations holding for the soft-overlap contributions are not 
broken by the hard-scattering terms, and therefore, at leading order in 1/E and to all orders 
in as, only eight of the ten form factors describing B P,V decays are independent. The 
relevant matching coefficients (jet functions) can be related to two universal quantities jTji 
and J'±, which we compute including one-loop radiative corrections. As an application of our 
results, we present in Section 7 a RG-improved analysis of the hard-scattering contributions 
to the large-recoil heavy-to-light form factors, in which all single and double logarithms arc 
resummcd. Section 8 treats the asymptotic limit of Sudakov suppression factors in the soft- 
overlap and hard-scattering terms. In Section 9 we present our summary and conclusions. 



2 Soft-collinear effective theory 

In processes involving energetic light particles, such as the pion emitted at large recoil in 
semileptonic B decay, it is convenient to introduce light-cone coordinates 

p'^ = n ■ p — + n ■ p — + p'^= p1 + p't + p'^, (2) 

the second equality serving to introduce the vectors p+ and p-. The light-like vectors n^, fi^ 
satisfy = = and n ■ n = 2. As mentioned in the Introduction, depending on the value 
of the rcnormalization scale the effective theory is described by hard-collinear and soft modes 
(SCETi), or by coUinear, soft, and soft-collinear messenger modes (SCETn). It is conventional 
to quote the scaling behavior of the components (p^,p^,p±) with the energy in terms of a small 
parameter A ~ Aqcd/^'- The coUinear and soft momenta of the partons inside the external 
meson states in 5 — > M transitions scale like p^ ^ E{X'^, 1, A) and p^ ~ E{X,\,\). Hard- 
collinear momenta arc defined to scale as phc ~ -^'(A, 1, A^/^), whereas soft-collinear momenta 
scale as Psc ~ E^X^, A, A'^^^). Throughout this paper we will identify 

E ^ V ■ p- ^^{n ■ v){n ■ p) Aqcb (3) 

with the energy of a coUinear or hard-collinear particle in the 5-mcson rest frame. Strictly 
speaking, the energy v ■ p differs from E by terms which can be neglected at leading and first 
subleading power. 

Fields in SCETj and their scalings with the expansion parameter A^/^ are [2, 3, 4] 

Al ~ (A, 1, AV2) , ~ (A, A, A) , (4) 

where h is the heavy-quark field in Heavy-Quark Effective Theory (HQET) [16]. The leading- 
order soft and hard-collinear quark Lagrangians are 

^hl = ^hc ^ [in ■ Dhc+s + iphc± .- iphc±] ^hc ■ (5) 
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Here iD^ — id^ + gA^^ denotes the covariant derivative built with soft gluon fields, etc. For 
the pure gauge sector we may write Cf^ = £5 + C^l^l , where Cs is the usual gluon Lagrangian 
restricted to soft fields, and is obtained by substituting iD^ — > ^-^L + 9^^+ ii^t*^ ^^c 
Yang-Mills Lagrangian. In interactions involving both soft and hard-collinear fields, the soft 
fields are evaluated at X-. In particular, this rule applies to the first term in the hard-collinear 
Lagrangian, in which in • Dhc+s — in ■ d + gn • Ahc{x) + gn ■ As{x-). Corrections from the 
multipole expansion (f)s{x) — (j)s{x-) +x±-d± (l)s{x-) + . . . appear as higher-order terms in the 
expansion in A^/^. We shall also need the subleading interaction Lagrangian [4] 

4f ^ = tphc±Whc Qs + h.C. , (6) 

which transforms a soft quark into a hard-collinear one (and vice versa). Here Whc is a hard- 
collinear Wilson line in the n direction necessary to ensure gauge invariance. Additional terms 
in the Lagrangian are required to describe the soft-overlap contribution in SCETj [6]; however, 
they will not be of relevance to our discussion here. 

The Lagrangian interactions are a special case of general gauge-invariant operators. Im- 
posing homogeneous gauge transformations in the soft and hard-collinear sectors [17], which 
strictly respect the SCET power counting, the "homogenized" hard-collinear fields are re- 
stricted to appear in the combinations 

X,e = W^iac~AV2, Al^^Wl{iDl^Wf^)r.xy\ wlm-D,,^,WHc-^X, (7) 

and may be acted on by partial derivatives in • d ^ 1 and id^ ~ X^^'^. This result follows 
from considering the most general operator in the gauge where n ■ A^c — 0, n ■ Ag — 0, and 
then returning to an arbitrary gauge by using the transformation laws appropriate for the 
homogenized fields. The 0{1) components in ■ d may appear an arbitrary number of times 
in operators of a given order in the power counting. This is accounted for by smearing hard- 
collinear fields along the n direction, i.e., (phcif^n) with arbitrary r. Soft fields may appear as 
Qs, h ^ A^/^ and iD^^ ~ A, and position arguments x±_ ~ A~^/^ and x+ ~ A° arising from the 
multipole expansion may also appear in interactions involving both soft and hard-collinear 
fields. Until Section 6 we deal exclusively with SCETj, and from now on will drop the label 
"/ic" on the hard-collinear fields. 



3 Flavor-changing currents in SCETj 

Heavy-to-light form factors describing current-induced S — > M transitions (with M a hght 
meson) at large recoil are subleading quantities in the large-energy limit, in the sense that the 
transitions they describe cannot be mediated by leading-order SCETj currents and Lagrangian 
interactions. For a leading-order analysis of these form factors, it is sufficient to include heavy- 
coUinear current operators through the first subleading order in SCETi power counting [6, 7]. 
These operators contain a heavy-quark field, a hard-collinear quark field, and (in the case of 
subleading operators) transverse derivatives and gluon fields. They provide a representation 
in SCETi of the QCD heavy-light current operators Jqcd = ^F 6, which we renormalize at 
a scale a*qcd- (The vector and axial-vector currents are not renormalized.) At a given order 
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in power counting, a minimal basis is determined by first writing the most general gauge- 
invariant operators constructed from the available fields and external parameters, and then 
requiring invariance under small variations of the external parameters [3, 18]. 

3.1 Determination of the operator basis 

It is convenient to restrict attention to the case = 0, where v denotes the S-meson velocity. 
It follows that 

n>' = ^(2v>'- . (8) 

n • V \ n ■ V J 

With this choice there are two remaining reparameterization transformations, which we con- 
sider in their infinitesimal form. The first enforces invariance under the rescaling 

n'' ^ (1 + a) n'^ , ^ {1 - a) (with a ~ 1) . (9) 

The second allows small changes in the perpendicular components, such that 

n^^n'^ + e''., ^ - , ^ (with - A^/^) . (lo) 

[n • vy 

The power counting assigned to a and is the largest possible (i.e., providing the strongest 
constraints) such that the scaling of hard-coUinear momenta is unaltered. In both cases, the 
variation of is determined by (8) and the variation of n^. An alternative approach would 
be to introduce n and v as arbitrary vectors, subject only to the conditions = 0, n ■ n = 2, 
and = 1. Then n'^ would not transform under (10), and there would be an additional 
reparameterization transformation 

n'^^n'', n" ^n^^ (with ~ 1) . (11) 

In this case a consistent power counting requires v±_ ~ C(l), and many more operators appear 
at a given order than when = [19] . At higher order it would also be necessary to impose 
invariance under heavy-quark velocity transformations, v ^ v + 5v with v ■ 5v = 0. However, 
these transformations enter only at C(A) and so are irrelevant for our leading-order analysis 
of heavy-to-light form factors. 

Requiring invariance under these transformations, it is straightforward to write down the 
most general operators with given quantum numbers. The leading-order currents are C(A^), 
and for the scalar case we find 

4°) (s, x) = e-^"^"'"^ X{x + sn) h{x.) , (12) 

where the phase factor arises from the definition of the HQET field h, and only the first term, 
h{x-), has been retained in the multipole expansion. In the form-factor analysis we can use 
translational invariance to set x = in the weak current operators, and so we restrict attention 
to Js^\s) = Jg^\s,x — 0). Similarly, for the vector and tensor currents we have 

4^^{s) = X{sn) Tfh{0) , 4'^^''{s) = X{sn) h{0) , (13) 
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with the Dirac structures 



n ■ V 



rf = 7[v', rf = wV', rf = ^--^, rr = — — . (i4) 

Square brackets around indices denote anti-symmetrization. Note that all structures invariant 
under the rescaling (9) are allowed at leading order, since the reparameterization transforma- 
tions (10) enter only at subleading order. 

Before writing the most general set of subleading operators, we first consider the variation 
of the leading operators under (10). To first order in A^/^, the fermion fields transform as 

SX{x + sn) = X{x + sn) - e_L • gAs±_{xJ)^ , 

5h{xJ) ^^e^- d^h{xJ) . (15) 

The soft gluon appearing in the variation of X(a;) arises from the field redefinition enforcing 
homogeneous gauge transformations in the hard-coUinear sector [17]. The combination X of 
hard-coUinear fields may be expressed as X = i?'^ (^.^74) ■j/'l^^, where is the (unhomogenized) 
hard-collinear fermion field in the gauge n ■ A^^ = 0, and R{x) is a soft Wilson line from x_ to 
X. Under the transformation (10), 5^^^ ~ ^(^)) the remaining variation of X arises from 
the projection ^^/4 and the homogenizing factor R. Next, we consider the variations of the 
various Dirac structures under the transformation (10), finding 



n ■ V 



5rf = o, 5rf = o, 5rf = ^^, sr^"^^^. (i6) 

n - V n - V 

The variations of the subleading operators must cancel these contributions. 

The subleading C(A^/^) operators must contain exactly one insertion of or {—id±^) 
acting on the hard-collinear field X, or x± • Ds± with the derivative acting on h. To determine 
the most general form, we note the following transformation properties, working now to zeroth 
order in A^^^: 

5A^^ = 0, s(^^)=-'^, s{x^-D,^) = -^e^.D,^. (17) 
\—in ■ d J ^ ^ 

Thus, (— i dD and x± ■ Ds± are restricted to appear in specific reparameterization-invariant 
combinations with the leading-order currents, and there are no constraints on the appearance 
of A'^. Inspection of (15), (16), and (17) allows us to deduce the form of the most general 
operators through 0(A^/^). For the scalar current 

Ji{s, x) = e-^"^*-^-^ X{x + sn) (l - 4^ - + x^- D,j^ h{x_) , 

V in - d J 
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J|(s, r, x) = 6-^"^"'"=^ X{x + sn) 4±{x + rn) h{x^) . (18) 

Again, using translational invariance we can specialize to x — and define Js{s) = J^{s,0) 
and Js{s, r) = J5 (s, r, 0). Similarly, for the vector and tensor currents at x = we obtain 



4r^,{s) = x{sn) (1 - M= I) r^,, MO) , 

V in - o / 



4,^{s) = X{sn) ( 1 - I ] /i(O) + ^ Xisn) h{0) , 



and 



J^^'{s, r) = X(sfl) yi^;;(rn) h{0) , (19) 

J^,'|'(.) = Xisn) (1 - ^ I) r^^ MO) , 

\ in ■ ^ / 

4r{s) = X(sn) fx - ^ ) rr /i(0) + ^ X(sri) /i(0) , 

\ in ■ d ^ J '^■'^ in ■ d 

4r{s) = X(sn) f 1 - ^\ rf /i(0) + ^ X(sn) /i(0) , 

V in ■ d ^ / ^■'^ in ■ d 

Jti,Iza{s, r) = X{sn) 4±{rn) r5'2,3,4 ^(0) , 

J?,7,7i', r) = X(srz) yit(rn) r^,,3 /i(0) . (20) 

In the tensor case, one combination of the J^i^'^ is redundant in four dimensions, being propor- 
tional to the anti-symmetric product 7|['7X7i'- This combination will be isolated in the next 
section, when we define a new basis of current operators that are renormalized multiplicatively. 

We will refer to the quark- ant iquark operators as "A-type currents", and to the quark- 
ant iquark-gluon operators as "-B-type currents". The A-type currents contain both leading 
and sublcading contributions, which are linked by reparameterization invariance. Pseudoscalar 
and axial- vector currents are obtained from the above expressions for scalar and vector currents 
by insertion of 75 next to the light fermion field. (This simple prescription holds only in the 
"naive dimensional regularization" (NDR) scheme.) The results (18), (19), and (20) agree 
with the basis of operators found in [19], specialized to the case where v± = 0. However, the 
derivation presented here is much simpler. We also stress that our definition of the subleading 
currents, in which the A-type currents contain d_i rather than a covariant derivative, will 
ensure that the A-type (two-particle) and S-type (three-particle) operators do not mix under 
r enormalization . 

The subleading contributions to the A-type currents containing perpendicular derivatives 
on the hard-coUinear fields do not contribute at leading order in the form factor analysis. 
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QCD SCETi SCETii 




Figure 1: Two-step matching procedure for a hard-scattering contribution. The dashed hues 

in the SCETj diagram represent hard-colhnear fields, those in the SCETn four-quark operator 
denote colhnear fields. In all cases, the heavy quark is depicted as a double line, and the wavy 
line represents the flavor- changing current. 



The extra derivatives yield an 0{X) suppression on top of the suppression factors already 
present when the leading currents mediate the decay process by either the soft-overlap or 
hard-scattering mechanisms. A formal demonstration of this point can be found in [7]. The 
remaining 74-type currents (at x — 0) all take the form X{sn)Th{0), differing only by their 
Dirac structure F. Spin symmetries arising from the constraints — 0, — may then be 
used to relate matrix elements involving the same initial- and final-state mesons. Although the 
SCETii representation of the A- type currents is rather complicated [8], owing to the symmetry 
relations the results can be expressed in terms of a small number of non-perturbative functions. 
The A- type currents thus give rise to the first, soft-overlap term in (1). 

In general, the B-type currents break the spin symmetries. When applied to form-factor 
matrix elements, the hard-colhnear gluon emitted from the current is absorbed by the spectator 
quark, allowing the decay to proceed via the hard-scattering mechanism. Figure 1 illustrates 
the two-step matching of a typical hard-scattering amplitude. In the first step, the QCD 
current is matched onto a B-type SCETj current. In the second step, the hard-colhnear 
gluon is integrated out, and the hard-scattering amplitude is described in the final low-energy 
theory by non-local four-quark operators. Section 6 examines this SCETn representation, 
where matrix elements take the form of the second, symmetry-breaking term in (1). 

3.2 Matching calculations 

We proceed to find the SCETi representations of the QCD scalar, vector, and tensor currents 

S^qb, V^q^^b, T^"" ^ {-i)qa'"'b^q-f^^'Y^b. (21) 

The QCD operators S and T'^'^ require renormalization and are defined in the modified minimal 
subtraction (MS) scheme at a fixed scale //qcd = 0{mb). The SCETj representations of these 
operators, evaluated at position x — 0, are given by an expansion (summed over i,j) 

qTb^ JdsCfis)jf{s) + ^ JdrdsCf{s,r)jf{s,r) + ... 

= Cf{E) Jf{0) + duCf{E,u) Jf{u) + . . . (22) 
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with operators Jf, Jf as determined in the previous section for the appropriate quantum 

numbers. We denote coefficient functions in position space with a tilde. In the second hne, 
we have used translational invariance and defined the momentum-space coefficients (without 
a tilde) as 

Cf{E,u) = J drdse'^'''+'''^''-^Cf{s,r). (23) 

Here P — Pout — -Pin is the total hard-coUinear momentum of external states (strictly speaking 
this is a momentum operator), and E — v ■ P- — {n ■ v){n • P)/2. Reparameterization 
invariance ensures that the momentum-space coefficient functions depend on the combination 
(n ■ v){n ■ P) = 2E, not n ■ P. The variable u G [0, 1] is the fraction of the large momentum 
component n ■ P carried by the fields in X (the dressed outgoing hard-coUinear quark field), 
and u = 1—uis the corresponding momentum fraction carried by the fields in A± (the dressed 
outgoing hard-coUinear gluon field). The object Jf{u) in (22) denotes the Fourier-transformed 
current operator 

— e-^"^-^^ Jf (s, 0) . (24) 

271 ■' 

Note that both Jf{Q) and Jf{u) also depend on the large energy scale E as well as on the 
renormalization scale /i. These dependences will be suppressed for simplicity. 

The matching conditions for the momentum-space Wilson coefficient functions Cf and 
at tree level follow from an analysis of current matrix elements in QCD and SCETj. The 
only subtlety is that a non-zero matching contribution is obtained from graphs where a hard- 
coUinear gluon is emitted from the hard-coUinear quark line. This might seem surprising at 
first sight, because the resulting propagator is close to the mass-shell. This contribution is 
present because two of the four components of the hard-coUinear quark spinor are removed 
when QCD is matched onto SCETj. To see how it arises, consider the following diagram: 



» ' — > — > = »' — » ts — > + ^ --» 

o! a>{ cpi 

o! cb< cK 

"=3 i i (25) 

7^:ir = — r + 7/.^^r 

p 2 n ■ p 2 n ■ p 

For vanishing transverse momentum {p± = 0), the intermediate quark propagator takes the 
form shown in the second line. In the effective theory, the first term on the right hand side is 
represented by a graph where the gluon is emitted from a hard-coUinear quark line. The second 
one, however, corresponds to a graph where the gluon is emitted from the B-type current, and 
contributes to the matching coefficient. Our tree-level results are given as follows. 

Scalar current: 

Ci = l, Ci = -1. (26) 
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Figure 2: One-loop QCD diagrams contributing to the matching calculation for the subleading 
scalar current. The external gluon can be attached at any of the places marked by a cross. 



Vector current: 




v\ — 



1, 



Cy2 — —2 , 





(27) 



Tensor current: 




a 



T2 — 



T2,3,4 



C't3 — ~2x , C'r4,5,6 — > ^T7 — ■ (28) 



Here x — 2E/mb. Our tree-level matching coefficients agree with previous results, for Cf in 



For the resummation of the leading logarithms, the tree-level Wilson coefficients are suffi- 
cient. However, for the physical value of the 6-quark mass the one- loop matching corrections 
may turn out to be comparable to the effect of leading-order running. To address this ques- 
tion, we evaluate as an example the one-loop matching corrections for the scalar current. We 
calculate the decay of a 6-quark to an energetic light quark and an energetic gluon at one-loop 
order. Throughout this paper we use dimensional regularization with d = 4 — 2e dimensions 
and employ the MS scheme to remove ultra-violet (UV) singularities. We use the background- 
field method and perform the calculation in an arbitrary covariant gauge. The eight one-loop 
diagrams contributing to the matrix element are shown in Figure 2. Each diagram involves 
at least one heavy-quark propagator; diagrams involving only massless propagators vanish in 
dimensional regularization once the external lines are put on the mass shell. Note that there is 
again a contribution in Figure 2 corresponding to gluon emission from a nearly on-shell hard- 
coUinear quark line, which is decomposed as in (25). To correctly identify the two parts of the 
QCD loop diagram, one first considers it for 7^ and then expands around = 0. This 
expansion must be done before performing loop integrations. This ensures that all effective- 
theory loop diagrams vanish and only the hard part of the QCD diagram is left, which is 
exactly the part that has to be absorbed into the Wilson coefficients. We set all perpendicular 
components of the external momenta to zero and equate the effective-theory expression to the 
QCD result for the three-point function. Depending on the polarization of the background 
gluon field, we thus determine either the Wilson coefficient Cg of the leading-order current 

^Notc that wo use square brackets around two or more indices to denote antisymmetrization (^j^f^j'^^ = 
^il'^l" ^ 1^1^), etc.), whereas in [19] square brackets around two indices denote a commutator (^I'^'y"! = 
7^7"^ -7'^7^, etc.). 



[1], and for C/ 



in [4, 19]. 
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operator (for an A_ gluon), or the coefficient Cg of tfie subleading current (for an A± gluon). 
After MS subtractions, we obtain in tfie first case (again witfi x — 2E/mb) 



Ci{E) = 1 + 



An 



-21n^ — + 51n— + 61n^ + ^ 



1^ 



2E 



— Inx — 2 Lioii — x) — — 
X ^ 12 



(29) 

where ag = c^sil-i')- The fact that this expression, which was extracted from a three-point 
function, agrees with the result obtained from the heavy-to-hght two-point function [1] follows 
from invariance under coUinear gauge transformations. Performing the calculation for the case 
of perpendicular gluon polarization, we obtain for the coefficient of the operator Jg 



C^{E,u) = -Ci{E) + 




4 + 



2x 



Inx 



X 



X 



1 — xu 



1 + 



1 



+ 



X 



1 — xu 1 — X 



\n{xu) 



+ 



2(2 -x) \nu 



1 — X 



u 




. , , 21nM 
l + \nu] — h 



2x 



u 



xu 



\n{xu) 



2\nu 



u 



+ 



xuu 



'1 — xu) 



Li2(l — x) — Li2(l - xu) — Li2(l - xu) + — 





71 



21 



(30) 



where u corresponds to the fraction of the longitudinal hard-collinear momentum carried by the 
final-state quark. We note that at the endpoints, Cg diverges only logarithmically: Cg ~ In^ u 
ai u — and Cg ~ IniZ at = 0. Also, despite appearances, there is no singularity at x = 1 or 
X — 1/u. The scale dependence of Cg agrees with a direct analysis of operator renormahzation 
given in the following section. We postpone a detailed discussion on the relative importance 
of matching and RG running until Section 5, when a solution to the RG equations is at hand. 

As a final remark before ending this section, we note that the scalar current is not an 
independent operator but can be related to the vector current using the equation of motion 
for the quark fields, idiV — fUbdiqcB) S (iiqcb) , where 5'(/xqcd) denotes the scalar QCD 
current in (21) renormalized at scale /xqcd, and m{,(/iQCD) is the running 6-quark mass, both 
defined in the MS scheme. Applying this identity to (22), it follows that to all orders in ag 
and at leading order in Aqcd/^6 



-]C^ 

rrib 



V2 



,r<A_ ^fe(/^QCD) 



nib 



1 - 



rribj 



V2 



+ Cv3 — 



mbifJ^QCD) 



rub 



s > 



(31) 



where at this order there is no difference between the meson mass ms and the 6-quark pole 
mass rub- It is readily seen from (26) and (27) that these relations hold at tree level. At 
one-loop order 



rrib 



1 + 



An 



6 In 



rUb 

/^QCD 



-4 



(32) 
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Note that the dependence on the scale Hqcd cancels between the running mass ^^(/xqcd) and 
the scalar coefficients Cg and Cg on the right-hand side of the relations in (31). The ffist of 
these relations can be verified to hold at one- loop order using the results of [1]. In Section 7, 
we will use these results to deduce the one-loop matching coefficients relevant to the vector 
form factor Fo(g^), using only scalar-current matching calculations. 

4 Anomalous dimensions 

The matching results derived in the previous section can be trusted at a high renormalization 
point ^ = 2E, at which the Wilson coefficients are free of large logarithms and so can 

be reliably computed using fixed-order perturbation theory. In order to evolve the coefficients 
down to lower values of /x, one needs to solve the RG equation for the SCETj current operators. 

The currents Jf and Jf do not mix under renormalization. To see this, we first note that 
the operators cannot mix into Jf, as they are of higher order in power counting than 
the leading terms of jf. In principle, the operators jf could mix into via time-ordered 
products of either the leading terms of Jf^ with the subleading SCETj Lagrangian, or the 
subleading terms of Jf with the leading SCETi Lagrangian. However, when v'^ — 0, in both 
cases the resulting time-ordered products have a structure different from the operators, 
containing additional perpendicular partial derivatives acting on hard-collincar fields. Finally, 
the various operators do not mix among themselves, since time-ordered products with the 
SCETi Lagrangian cannot mix the different Dirac structures appearing in their definition. As 
a result, the effective-theory currents obey the integro-differential RG equations (summation 
over k is understood in the second equation) 

7^J^(0), 

fdv^%MJ^{v), (33) 
and the corresponding momentum-space Wilson coefficients satisfy the equations 

fdu^ij{u,v)Ci{E,u). (34) 

The anomalous dimensions 7^ and 7^^ are calculated by isolating the UV divergences in SCETi 
loop diagrams. The one-loop expression for 7"^ has been calculated previously [1], with the 
result 

7^ = -rc.sp(«.) In ^ + 7(a.) = ^ (" 1^ ^ " + ^(«') ■ (35) 

In [20] it was argued that the first equality in (35) is valid to all orders in a^. The appearance 
of In n in the anomalous dimension is explained by the theory of light-like Wilson loops with 



d 



dln/j, 
d 

din a 



^^(0) = - 
Jfiu) = - 



d 



din// 

d 



cm 



dmn ■' 
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cusps. Only a single logarithm appears at any order in the strong coupling, with a coefficient 

governed by the universal cusp anomalous dimension Fcusp [21]. 

Unlike the situation for the A-tjpc operators, the anomalous dimensions of the B-tjpe 
currents depend on the specific Dirac structure, and there is non-trivial operator mixing. The 
renormalization properties of all B-type operators can be discussed by defining two operators 



(36) 



and their Fourier transforms as defined in (24). The Dirac structure is determined by the 
specific current under consideration. It follows from the Feynman rules of SCETi and the 
projection properties of the two-component spinor X that the two operators 3i close under 
renormahzation. Defining renormahzation constants via 3]^"^ — Z^j d^'^^^, we obtain 



Zn Zi2 
Zu + 2{l-e)Zu 



(37) 



poles 



where, in the MS scheme, it is understood that only pole terms in the dimensional regulator 
e — 2 — d/2 are kept in Z — 1. We define two anomalous-dimension functions 



7i = 2a 



d 



7{l) 
'11 ) 



72 = 2as 



d 



.(1) 



' das ' " " das ' 



(38) 



where Z^^^ is the coefficient of the 1/e pole in Z. Then, at one-loop order, the corresponding 
anomalous dimension matrix determining the running of the currents 3i reads 



'Vl— loop 




(39) 



At higher order, the 0{e) terms arising from the Dirac algebra have a non-trivial effect on the 
anomalous dimensions, as discussed in [22]. 

In terms of these functions, the anomalous dimension of the scalar current is (from now 
on all results refer to the one-loop approximation) 



7f 



7i + 272 ■ 

For the vector current, we find the 4x4 anomalous-dimension matrix 



(40) 



71 + 272 
71 + 272 
71 + 272 



\ 72 







-72 








(41) 



71 / 



The operators Jyi2 3 ^^'^ multiplicatively renormalized with anomalous dimension 71 + 272, 
whereas the operator Jy^ mixes with Jy-^ and Jy^. For the form factor analysis, it will be 
convenient to work with the linear combinations 



[rn 



ilKo), 
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jB'/j, 


_ jB/j. _ 


X(, 






X(, 


jB'ij. 


_ 9 





sn) A±(rn) h(0) , 

n ■ V 

+ J?s' = 4±irn) h{0) , (42) 

which are multiphcatively renormahzed. Jyj 2,3 have anomalous dimension 71 + 272, while Jy^ 
has anomalous dimension 71. The corresponding combinations of Wilson coefficients, which are 
multiphcatively renormahzed with the same anomalous dimensions as the respective currents 
read 

1 
2 



^B' _ ^B I ^ /^B r-iB' _ r-iB r^B' _ r^B , r^B 



C^. = \C^.- (43) 

Similarly, for the tensor current we find that the operators J^^ 234 multiphcatively renor- 
mahzed with anomalous dimension 71 + 272, whereas the operators Jtq^j with J-^i 2,3,4- 
As in the vector case, it will be convenient to change the operator basis, defining 

J^;'^'^ = J^f'^ + 24 f - = "^{sn) 4^{rn) 7^ ^(0) , 

J^'^^'- = J^/- + J^/- = X(sn) 4^{rn) v^^^'l h{0) , 

'^Ts'^'^ = Jx3^'^ = X(sn) ^j_(rn) — ^ /i(0) , 

Jt/'" = <^T4'^'' = X(sn) v^j_(rn) /i(0) , 

n ■ V 

J^V = -24^^" + 41"" + 2 J^g'^'^ - 2 J^/'^ + 2 J^,'^'' = l{sn) AUrn) h{Q) , 

J^^''- = -2 J^e'^" - J^/" - J^/" = X(sn) ^;1'^7^ ^^(m) /i(0) , 

[a* 

J^^'^'^ _ _2 j^^/^- _ 7^3^^- = X(sn) ^^-^ ^x(rn) /i(0) . (44) 

77/ * V 

These operators are multiphcatively renormahzed, ^^1,2,3,4 with anomalous dimension 71 + 
272, and Jt5,6,7 with anomalous dimension 71. The corresponding combinations of Wilson 
coefficients are 

^^B' _ ^B I 1 ^B ^B' _ ^B 1 ^B /~iB' _ ^B on^ ^ ^ 

^Tl — ^Tl "T 2 1 ^T2 — ^T2 ^ '^6 ' ^TZ — ^T?, ^^Tl ^ '^^ 2 ' 



^T4 — T4 ^'^Tl T2 ' 
^B' _ 1 ^B r~iB' _ _ 1 r~iB ^B' _ _ 1 _ 1 //I f-\ 

Note that the tree-level matrix elements of the "evanescent" operator J^^ vanish in c? = 4 
dimensions, and after additional finite renormalizations the same is true once loop corrections 
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Figure 3: SCETj graphs contributing to the anomalous dimensions of the subleading heavy- 
coUinear currents , represented by a crossed circle. Pull lines denote soft fields, dashed lines 
hard-coUinear fields. 

are included. This operator will not be relevant for our leading-order analysis, but it would 
enter a next-to-leading order calculation employing dimensional regularization. We will men- 
tion evanescent operators again in Section 6, when we discuss the basis of SCETn four-quark 
operators relevant to the hard-scattering contributions. 

To calculate the one-loop anomalous dimensions, we evaluate the UV poles of the diagrams 
shown in Figure 3 in dimensional regularization, treating the external gluon as a background 
field. The resulting expressions for the integral operators 71 and 72 can be written as 



ji{u, v) — u Vi{u, v) -\- 5{u — v) W{E, u) , 
272 (1*,^;) = uV2{u,v) , 



(46) 



where 



Vi{u,v) 



V2{U,V) 



a. 



2 / TT 

1 



d{l-u-v) 



UV 



Ca Pis 

2 t: \ UV 

2 TT 



{u — v) 
{v — u) ' {u — v) 



e(v - u) 

U^ ^ + V' 



9{v — u) 9{u — v) 



1 1 



vu 



uv 



UV 



i + - + -\ eii-u-v) + — eiu + v-i) 

U V } UV 



Ca_OLs_ 
2 TT 



I + V\ eiu - v) + - (1 + \ ] eiv - u) 

V V \ u ' 



W{E,u) = C,^ (-X^^-U\nu\ - fln^ - l) . 

^'^ TtV 2E a J 2 TT \ u J 



(47) 



For symmetric functions g{u,v) the plus distribution is defined to act on test functions f{v) 



as 



dv[g{u,v)]+f{v) = / dvg{u,v) [f{v)-f{u)] 



(4J 
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It is remarkable that the functions Vi{u, v) are symmetric in their arguments. This fact is not 
accidental, but can be traced back to a residual conformal symmetry in the effective theory 
when interactions with the soft heavy quark are ignored. Also, since soft gluons are unable 
to change the large component of hard-coUinear momenta, the one-loop conformal-symmetry 
breaking term 5{u — v)W{E, u) is local with respect to u. Details of the conformal-symmetry 
arguments will be presented in Section 5.3. 

From (33), it follows that the evolution equation for the current eigenvectors jf takes the 
form 

d 



Jf\u) - - / dvuVr{u,v) Jf\v) - W{E,u) jf\u) , (49) 



where is a linear combination of Vi and V2 determined by the anomalous dimension of the 
operator . For reasons that will become clear later, we use the notation F = ||, _L and denote 

V|| = + V2 for operators with anomalous dimension 71 + 272, and Vj_ = Vi for operators 
with anomalous dimension 71. The corresponding equation for the eigenvectors of Wilson 
coefficients reads 



d 



^, Cf\E,v)^ I duuVr{u,v)Cf{E,u) + W{E,v)Cf{E,v). (50) 
dhiii ' Jq ^ 

To gain more insight into the structure of the conformal-symmetry breaking term W{E^u)^ 
we may consider the field redefinitions [2, 15] 

inlx) = S^{x.) el?(a;) , Al{x) = S^{x.) Af/{x) SUx.) , h{x) = S^{x) h^'\x) , 

(51) 

which decouple soft gluons from the leading-order hard-coUinear and heavy-quark Lagrangians. 
Here Sn and are soft Wilson lines in the n and v directions, respectively. From (7), it follows 
that in terms of the new fields the i?-type operators take the form 

Jf(s, r) = X^'\sn) Aflirn) T'i [Sl S,] (0) h^'\0) . (52) 

The combination [5*^ Sy] (0) represents a closed loop with a cusp at x = formed by two 
Wilson fines in the v and n directions. The anomalous dimension of this object is given by 
the universal cusp anomalous dimension times a logarithm of the soft scale [21]. After adding 
a contribution from the hard-collinear sector necessary to eliminate the dependence on the 
infra-red regulator {—pic) (c-f- [15, 23],), the result is 



r 

^ cusp 



■^^2t^^^^(-pL) 



i-plc) ^^ 



2 



Tcusp In , (53) 

A* 



with a hard-collinear momentum. From these considerations, we conclude that to all orders 
in perturbation theory 

W{E,u) = -rcusp(Qi.) In ^ + w{u, a,), (54) 
with the one- loop expression for w determined from (47). 
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5 Renormalizat ion- group evolution in SCETi 



The anomalous dimensions obtained in the previous section allow us to solve the RG evolution 
equations (49) and (50) at leading order in RG-improved perturbation theory. At this order the 
leading double and single logarithmic contributions are resummed to all orders in perturbation 
theory. Technically, this means that one must compute the matching coefficients at tree level, 
the anomalous dimension kernels Vi{u, v) at one-loop order, and the cusp anomalous dimension 
entering the function W{E,u) at two-loop order. In practice, one wants to impose matching 
conditions for the coefficient functions Cf at a high scale jih ~ '^E and then evolve the 
result down to an intermediate scale jii ~ y/2EKqcT>-, at which SCETi is matched onto the 
low-energy effective theory SCETh. While the integro-differential evolution equations can be 
solved numerically, we find it instructive to also discuss a formal analytic solution to these 
equations. 

5.1 Eigenf unctions and eigenvalues of the evolution operator 

It will be convenient to rewrite (49) and (50) in the somewhat obscure form 

— — / dvvv — 33 — - W{E,u) 



— ' dvvv — 33 — -— — z \-W{E,u)— — z ■ (55) 



dlnfi u Jq uv V u 

Because of the symmetry property Vi{u,v) = Vi{v,u), it follows that the operator eigen- 
functions [J^ {u)/uu] have the same form as the coefficient eigenfunctions [C^ {u)/u], but 
come with eigenvalues of the opposite sign. We consider first the hypothetical case where 
W{E,u) — and focus on the non-diagonal terms in the evolution kernels, Vi{u,v) with 
i — 1,2. We will show that the general solution to the eigenvalue equation 

/ du uv? ^^^'^l^^ tpn{u) = Aj „ ipn{v) (56) 
7o uv 

is given by 

Mu) = ^^^^^^pL»\2u-1), (57) 

where Pn ' are Jacobi (hyper-geometric) polynomials, and the eigenfunctions ipn are normal- 
ized according to 

{ipnl'lpm) = / duuu'^1pri{u)lprn{u) = 6nm- (58) 







The eigenvalues for the kernels Vi and V2 may be expressed in the closed form 



2 J TT n + 2 2 n \ ^ n + 2 
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V ^ 2 y TT (n + l)(n + 2)(n + 3) 2 tt (n + l)(n + 3) ' ^ ^ 

where — Y^m=i^/''^ harmonic numbers. Using the solution to the eigenvalue 

equation for Vi, we then present a formal algebraic solution to the evolution equation for the 
SCETi coefficient functions. 

Let us now prove the statements just made. We proceed along lines similar to the diag- 
onalization of the evolution kernel for the pion LCDA, as analyzed in [24] (see in particular 
Appendix D). Prom the eigenvalue equation (56) and the symmetry of Vi{u,v) it follows that 
eigenfunctions belonging to different eigenvalues must be orthogonal in the measure uu^ du on 
the interval < m < 1, as shown in (58). To show that the eigenfunctions take the form (57), 
we consider the integrals 



Jo 



We will show below that Ii n{v) is a polynomial in {2v — 1) of degree n. Hence, in operator 
notation, with {u\n) = {2u — !)"■ and (M|Vi|n) = Xj_„(ti), we have 

n 

Vi\n) = 5^(V^0mn|m) (61) 

m=0 

with {Vi)mn = ^ ioY m > n. In particular, |0) is an eigenf unction with eigenvalue (Vi)oo- By 
induction, for each n > 1 there is a corresponding eigenvalue given by the diagonal matrix 
clement (l^)nn with eigcnfunction proportional to the linear combination of |0), . . . ,\n) that 
is orthogonal to each of |0), . . . , |n — 1). The result (57) follows since the Jacobi polynomials 

(2 1) 

Pn ' {2u — 1) form the unique extension of the constant function to a basis of orthogonal 
polynomials with measure uu"^ du on the unit interval. As a byproduct of this analysis, the 
coefficient of {2v — 1)" in the expansion of 2i,n(v) is identified with the n-th eigenvalue Aj,„. 

To evaluate 2j^„, it is convenient to introduce new variables x — 2u — 1 and y — 2v — 1. 
The resulting integrals yield 

Ca\ 



= - ^) ^ (-1)-^ 



-Jn + l)(n + 2) 



2 TT 



n-1 



m=0 ^ ' ^ ' m=0 ^ ' ^ ' 



m 



C 1 " 

- f i 2(n + l)(n + 2)(n + 3) ^„ + ^"X-^)""" + ' + 2™] . (62) 



Inspection of these results shows that the coefficient of y"* in indeed vanishes for m > n, 
while the coefficient of gives the eigenvalues (59). 
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5.2 Solution to the evolution equation 

Due to the presence of the conformal-symmetry breaking term W{E,u), the eigenfunctions 
of the full anomalous-dimension operators 71 + 272 and 71 cannot be written in closed form. 
We will now show how a formal algebraic solution can be obtained. We begin by isolating 
the strong /i dependence of the coefficient functions due to Fpusp, writing the solution to the 
evolution equation (55) as 

Cf'{E,u,i^) = (—] [\yUr{u,v,i^n,l^)Cf' {E,v,iin) , (63) 

\l^h y Jo 

where 

with P{as) = das/dln/i, is a universal Sudakov factor, and 

/""'^^^^ da 

a(l^i,l^2)= rcusp(Q;) ■ (65) 

Note that S{fii,fi2) is negative, since rcusp(Q;s) > [25]. At leading order in RG-improved 
perturbation theory one finds [20] 



sii^h.iA 



An 



To 



a{l^h,l^) = -TT^ ^T^n, (66) 

where ri = as{f^)/as{fJ^h) > 1- The relevant RG coefficients arising in the perturbative expan- 
sion of the cusp anomalous dimension and the /3 function, 

°° +1 °° +1 

P{as)^-2asY,Pn{^y , r,„,p(a,)=5]r„(^)" , (67) 

n=0 n=0 



are To = 4C^, Ti = 4C^^ [(f - ^)C^ - f 7>ny], and /3o = ^ - fT^^n/, /3i = f Cl - 
"^CATpnj — ACpTpUf. We take n/ = 4 as the number of hght quark flavors, even at low 
renormalization scales. 

The remaining evolution is governed by a function Ur{u,v, fih, fJ') with initial condition 
Ur{u,v, Hh, IJ'h) = S{u — v) at the high-energy matching scale n = Hh- From (55) we find the 
corresponding RG equation 

— z = / dyyy — — z hw{u) = , 68) 

din n u Jq yu y u 

where 'w{u) is defined via (54). As before, the subscript F = ||, ± serves as a reminder that we 
must distinguish two cases of evolution functions corresponding to the different eigenvalues of 
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the anomalous-dimension matrices. A formal solution of this equation may be constructed by 
expanding Uy{u,v, /ih, l-i') on the basis of eigenfunctions ipni^), 



= ^ / ^ '^riW Un{v, IJ,h, 1^) , 



(69) 



where 



Un{v,IJ'h,IJ') = / 

Jo 



Ur{u,v,iJ,h,iJ,) 



u 



(70) 



follows from (58). Inserting this expansion into the evolution equation, and projecting out the 
coefficients Un, we obtain 



d 



din// 



(71) 



where Ar,n = ^i,n + ^2,n or Ar,n = ^i,n depending on the anomalous-dimension eigenvalue, and 
the quantities Wmn are given by the overlap integrals 



dvVv'^ljjm{v)w{v)llJn{v) ■ 



(72) 



Collecting the elements Wmn into an infinite-dimensional matrix w, and the eigenvalues Ar,i 
into the diagonal matrix Ar, we can write the formal solution to (71) in the form 



Pexp 



da 



[Ar(Q;) -\-w{a 



(73) 



where the initial condition at the matching scale /i = follows from (70).^ The symbol "P" 
denotes coupling- const ant ordering, with Ar(a;) + i«;(«) appearing to the left of Ar(tt') +'w{a') 
when a > a'. Combining this result with (69) provides a complete solution for the evolution 
function. 

For a leading-order solution, we may expand the anomalous dimensions and /3 function to 
one-loop order, defining as usual 



\ r \ \{o) , 



47r 



(74) 



Then the ordered exponential in (73) may be written as 



Pexp 



da 



as{lJ.h) 



P{a) 



[Xr{a)+w{a)] j — U exp 



Aiy^+t/)W)diag 

2/3o 



In ri 



U 



-1 



(75) 



^The expansion of {u)lu on the truncated basis of eigenfunctions {VAiCm), < n < N} converges in the 
hmit N ^ oo, with the norm (58), provided that j^duuu^ (u)/up < oo. Since it is only the convolution 
over the product of with the jet functions and meson LCDAs which appears in the final expression for the 
form factors, cf. (143) and (147), this condition is more restrictive than necessary on the endpoint behavior of 
the coefhcient functions, and this "norm" sense of convergence is stronger than we require. 
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Figure 4: Functions Ur{u, fi^, fi) in (77) for fi^ = 4.8 GeV and fi = 1.55 GeV. The curves 
correspond to RG evolution with 71 + 272 (dotted) and 71 (dashed). The left and right plots 
show results derived by using 20 and 40 basis functions ■0n, respectively. The solid hues were 
obtained by numerical integration of the evolution equation. 

where U is the unitary matrix that diagonahzes Ap-* + i.e. + to^^^ = C7 (Ap ^ + 

iy^*'^)diag Since A[?'* + to*^"^ is a real symmetric matrix it can be diagonalized with 

real eigenvalues, which we collect in the diagonal matrix (Ap^"* + if '■°-')diag- Finally, we can 
simplify the answer further by using that at tree level the initial conditions for the Wil- 
son coefficients at the matching scale /i = Hh collected in (26)-(28) are independent of u: 
Cf {E,u, jih) = Cf {E, jih). We then obtain the final result (vahd at leading order in RG- 
improved perturbation theory) 

Cf\E,u,^)^{^-£j^''^''^ e'Mu^^u,^n,^)Cf\E,^n), (76) 

where 

Ur{u, iih, IJ') = / dvUr{u,v,iih,IJ') (77) 
Jo 



= S V(" + l)(" + 2)(« + 3) 

In practice we will truncate the basis of eigenf unctions, so that Ap ^ + w^^^ becomes a finite- 
dimensional matrix, which can be diagonalized without much difficulty. 

Figure 4 illustrates the effects of RG evolution. We show numerical results for the two 
evolution functions Ur{u, Hh, fj) corresponding to the cases with anomalous dimensions 71 -|- 272 
(dotted curves) and 71 (dashed curves). The plots are obtained with = 2E — rUb — 4.8 GeV 
and ji — jii — \/2Eh.h — 1.55 GeV, where ^ 0.5 GeV serves as a typical hadronic scale. In 
addition to the solution obtained with 20 and 40 basis functions, the figure also shows results 
derived from a numerical integration of the evolution equation (solid lines). For the numerical 
solution, the evolution to a lower scale is performed in discrete steps of A In = 0.02. To 



t/exp 



.(0) 



+ w 



(oy 



diag 



2Po 



Inri 



U 
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calculate the change in Ur{u, /ih, /J'n) in the evolution step from /x„ to the convolution 

integral with the evolution kernel is evaluated for one hundred different u values, and the 
function Ur{u, fi^, fin^i) is obtained from a fit to these values. The results from the two 
different methods agree nicely: the curves obtained with a finite number of basis polynomials 
oscillate about the numerical results, the number of turning points being equal to the order 
of the highest basis polynomial. The amphtude of the oscillations decreases once more basis 
polynomials are included. To obtain the Wilson coefficients at the low scale, these results 
must still be multiplied with the universal Sudakov factor [2E / HhY'^^'^'^^ e^^^'^'^'' ~ 0.89 (for 
our choice of parameters). We observe that the additional, w-dependent resummation effects 
described by the functions Ur{u, /ih, Ij) are very small for the coefficients with anomalous 
dimension 71 + 272, whereas they are more sizeable for those with anomalous dimension 71, 
reaching about —20% for the smallest values of u. 

To summarize our results, we compile the Wilson coefficients for the -B-type current op- 
erators obtained at leading order in RG-improved perturbation theory. We introduce the 
short-hand notation Cr{E,u,n) = {2E/ |^hT'^^'^^'^^ e^^^'h'^'^ Ur{u, jih: 1^): where T =|| or T = ± 
depending on whether the anomalous dimension is 71 -|- 272 or 71, respectively. We will see later 
that these two cases are in one-to-one correspondence with the nature of the light final-state 
meson in B ^ P,V transitions. The case F = || applies for transitions into a pseudoscalar or 
longitudinally polarized vector meson (P or Vj|), while the case F = _L applies for transitions 
into a transversely polarized vector meson {V±). Our results are given as follows. 



Scalar current: 



Vector current: 



C7|(//) = -qi(£;,K,//). (78) 



CU^^) = -—C\\{E,u,^i), C^4(/^)=0. (79) 
rrib 

Tensor current: 

= -Cii{E,u,fx) , C^M = -4C||(E,«,/.) , C^4,5,6(/^) = 0, 

AE RE 

C^sif^) = C^{E, u, /.) , C^,{fi) = C^{E, u, /x) . (80) 

rub rub 

The corresponding Wilson coefficients in the primed basis follow from (43) and (45). At 
leading order, the results for the primed coefficients Cyi 2 4 and C^i 256 are given by the same 
expressions as in the original basis, i.e. = Cf, while Cy^{n) = (1 — 2E/mb) C\\{E,u, ji), 
and C^'^iiS) = C#4(/i) = 2C||(E,'u,/x), C^7(/x) = {AE/nib) C±{E,u, ji). 

Figure 5 illustrates our results for the scalar current. While one should use tree-level 
matching combined with one-loop running for a consistent treatment at leading order in RG- 
improved perturbation theory, the figure also shows results obtained by including the one- 
loop matching corrections presented in (30). We solve the evolution equation (50) by direct 
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Figure 5: Results for the Wilson coefficient Cg {E,u, /Ji) at E = mf,/2. The dashed Unes 
represent the tree-level (gray) and one-loop (black) coefficients at the high scale ii^ = 4.8 GeV. 
The sohd hnes are obtained after evolving both coefficients with the leading-order anomalous 
dimension down to the intermediate scale /li — 1.55 GeV. 



numerical integration, expanding the cusp anomalous dimension entering W{E, u) to two-loop 
order, but using one-loop expressions for Vi{u,v) and the remaining terms in W{E,u). We 
again take Hh = '2E = mi, and /Xj = ^lEK^ (with = 0.5 GeV). The dashed curves give 
the one-loop matching results at the high scale fih, while the solid lines show the result of 
RG evolution to the intermediate scale /Xj. Comparing the black and gray curves, we observe 
that the one-loop matching corrections are of the same order of magnitude as the effects of 
RG evolution. This fact provides motivation for an extension of our anomalous- dimension 
calculation to the two-loop order, which would be necessary for a systematic treatment at 
next-to-leading order in RG-improved perturbation theory. 

For completeness, we briefly discuss also the solution to the RG equation for the coefficients 
of the A-type currents, given in the first hue in (34). The solution is 



Cf{E,^l)=(—) exp 



. /' ~ Ct{E,^,n). (81) 



with S{iih,iJ,) and a{fj,h, fj) as defined in (64) and (65). At leading order in RG-improved 
perturbation theory we may use the expansions (66) together with 

n'-^^ da , 7o , 
P(a) 2po 

where 7(^5) is expanded similarly to (67), and 70 = —5Cf- As an illustration of the size 
of one-loop matching corrections, we may consider again the scalar case, where at tree level 
Cg = 1. At leading order in RG-improved perturbation theory, with tree-level matching, we 
find C^{E, jjii) Ri 1.097 at Hi = y/WK^ and E = mi,/ 2. From (29), the coefficient at the high 
scale = 2i? = m^, including one-loop corrections is Cg{E,fih) fa 0.934. Leading-order RG 
evolution to the intermediate scale then yields Cg{E,ij,i) a; 1.025. 
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5.3 Constraints from conformal symmetry 

The conformal invariance of QCD at the classical level can be used to simplify the solutions 
to the evolution equations for the hard-scattering kernels and hght-meson LCDAs. As a 
consequence of this approximate symmetry, the evolution equations become diagonal at leading 
order once they are written in a basis of eigenfunctions of definite conformal spin. In the case of 
the pion LCDA, these basis functions are the Gegenbauer polynomials C^"^ [24]. For the heavy- 
light current the situation is more complicated: the conformal symmetry is not only violated 
by quantum effects, but broken explicitly by the presence of the heavy quark. Interactions with 
the soft sector of the effective theory destroy the conformal invariance of the hard-collinear 
sector. However, since the soft interactions do not change the large momentum fractions of 
the hard-collinear particles, the breaking of the symmetry can only occur in the local part of 
the anomalous dimensions, i.e., in the term b{u — v) W{E^u) in (46). We focus here on the 
hard-collinear sector of the effective theory, showing in particular that the eigenfunctions of 
the non-local terms v) in (46) must take the form (57). For the remainder of this section 
we will refer to hard-collinear modes simply as "collinear" . 

Before discussing the properties of heavy-light currents in more detail, we briefly recall some 
aspects of conformal symmetry [26]. The full conformal algebra consists of the generators P^, 
3 of translations and Lorentz transformations, augmented by the generators D of scale 
transformations, — > Ax'*, and of special conformal transformations, 

^ x^ + a^x^ 

^ l + 2a-x + a2x2- ^^^^ 

The action of the generators on a field of scaling dimension I and arbitrary spin is 

i[P^, $(x)] = d^^{x) , i[J^„ ^{x)\ = (x^d, - x^d^ - E^,) ^x) , (84) 

i[D, ^x)] = {x-d + l) ^x) , t[K^, $(a;)] = {2x^x ■ d - x^d^ + 2lx^ - 2x^Y.^,) ^{x) . 
The spin operator E^,^ on scalar, fermion, and vector fields is given by 

S/iV = 0, E^i, 1p = l^^lJtu 4' 1 ^fiu = Qua ^^j, ~ 9fia ■ (85) 

The collinear part of the SCET operators is given by products of fields smeared along the light 
ray x'^ = rn'^, i.e. $(a;) = $(rn) = $(r). The SL(2, R) subgroup that maps this light-ray onto 
itself is called the collinear conformal group. It is obtained from the four generators n ■ P, 
n ■ K, n'^n^Ji^i,, and D. To classify operators under this group, it is convenient to work with 
the linear combinations 

L^i — —in ■ P , L_i = -n-K, 

Lo - '-(D + ^n^fC'J,,) , E='-{D- ^n^n"J,,) . (86) 

The generator E counts the twist t and commutes with the remaining three generators, which 
fulfill the angular-momentum commutation relations 

[L+,,L_,]=2Lo, [Lo,L±i] = ±L±^. (87) 
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The action of these operators on the fields $(r) is 



[X+i,$(r)] = -dr^r), 
[X-i, $(r)] = (r^ + 2jr) $(r) , 
[i:o,$(r)] = (r9,+i)$(r), (88) 

where the quantum number j is referred to as conformal spin. The twist t and the conformal 
spin j of a field $ are related to the scaling dimension I and the spin projection s through the 
relations t — I — s and j — ^{l + s), where s is defined by 

E$ = ^n'^n'^E^^^ = (89) 

To express the SCET fields in terms of (primary) conformal fields of definite spin and twist, 
we introduce the field-strength tensor 

gG^. = W^gGi^^W = i[D^, V,] , (90) 

where = — iA^. In particular, 

n • = n''g''/gSa/3 = gGL, n-dn-A^ n'^n^gQap = gQnn ■ (91) 

Then for the coUinear SCET fields at the origin, we find for the twist and conformal spin 
eigenvalues 

[£;,X(0)] = X(0), [Xo, X(0)] = X(0) , 

[E, g^Q)] = 2 g^Q) , [io, gnnm = ^;n„(0) . (92) 

In the following, we want to decompose a given operator into components with definite 
conformal spin. The construction of the corresponding basis is done in two steps. First, 
one identifies the operator of minimal conformal spin Oq^ i.e., the operator at the bottom of 
an irreducible conformal tower, defined by [L_i,(9o] = 0. The complete basis of conformal 
operators is then given by repeated application of the raising operator L^i to the highest- 
weight operator Oq.^ 

Ofc = [L+i,[X+i,[...,[X+i,Oo]...]]]. (93) 

^ V ' 

k times 

A trivial example of this procedure is the expansion of X(r) into operators with definite 
conformal spin. In this case the highest- weight operator is Oq = X(0), and the raising operator 
L^i acts as a derivative with respect to r. The expansion of X(r) in conformal spin thus 

^It is standard terminology to refer to the operator of minimal conformal spin as the highest- weight operator. 
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coincides with the Taylor expansion about r = 0. In order to analyze the SCET currents 
J^, we now consider the decomposition of the product of two fields, which we assume have 
individual conformal spins ji and j2, into operators of definite conformal spin. We first rewrite 
each of the two component fields <l>ji(ri) and $j2('"2) in the conformal spin basis. In other 
words, we expand the product in a Taylor series about ri — r2 — 0. At the n-th order, this 
leaves us with operators of the form 

5;^5;^%i(n)$,2(r2)|,^^,^^,, (94) 

where rii + n2 = n. In general, these operators do not have definite conformal spin. The 
minimal conformal spin of an operator built from two fields of conformal spin ji and j2 with n 
derivatives on the fields is j — ji + j2 + n, and the highest- weight operator for this case is [27] 




The Jacobi polynomials Pn appear as the Clebsch-Gordan coefficients of the coUinear con- 
formal group. 

As an illustration, the operators relevant for leading-twist light meson LCDAs have the 
form X(ri) r (^/2) X(r2). Prom the general result (95), and the conformal spin eigenvalues 
(92), it follows that the highest- weight operators are 

0("'i'i)(r) = {in ■ 9)"|x(m) T ^ P^^f " ^ ~ '"^ " 5 ^ X(rn)| . (96) 

I 2 \^^^ ■ d + in ■ d J ) 

For example, with F = 75, the pion-to- vacuum matrix element yields a moment of the pion 
LCDA, 

(0|O("'^'^)(0)|7r) oc fduPt'\u - u) Mu,^i) , (97) 
Jo 

which projects out the component of 0^(m,/x) proportional to uu Pri^'^\2u — 1), where P-i^'^^ oc 
Cf/'^- Similarly, we may expand the collinear fields yi^(ri)X(r2) appearing in the current 
operators into operators of definite conformal spin. For this case we find that the highest- 
weight operators are 

I \in ■ d + in ■ d J J 

and the corresponding eigenfunctions for are proportional to uu Pj^'^\2u — 1) , in agreement 
with (57). From the discussion following (55), the eigenfunctions for the Wilson coefficients 
must then be proportional to uPn^'^\2u — 1). 

In the conformal-symmetry limit L"^, Lq, and E are conserved charges, and operators with 
different conformal spin or twist therefore do not mix. For the collinear fields appearing in 
the heavy-light currents the situation is especially simple, since there is only one operator for 



26 



a given set of quantum numbers. Operators with different coUinear field content or witfi 9^ 

or n ■ d derivatives on one of the two fields have higher twist. Also, by construction, each one 
of the operators O^''^'^^ has unique quantum numbers, namely 

1^0,0^'"]= + (99) 

To the extent that conformal symmetry is preserved, the operators O^^!^''^'^^ are thus multi- 
plicatively renormalized. This explains why the eigenf unctions of the non-diagonal part of the 
leading-order evolution kernel take the form of Jacobi polynomials Pn , as was shown by 
explicit computation in Section 5.1. 



6 Renormalization-group evolution in SCETh 

After RG evolution down to the scale fii ~ ^y2EAQCD, the intermediate effective theory 
SCETi is matched onto SCETh. In this section, we present a complete operator basis for the 
hard-scattering contributions in SCETn and compute the corresponding matching coefficients 
(called jet functions). To investigate the size of loop corrections, we present the complete 
one-loop matching corrections to these coefficients and show that they can be expressed in 
terms of only two universal functions. We then consider resummation in SCETn and present 
the general solution to the RG equation. In the following Section 7, we will apply these results 
to obtain RG-improved expressions for the hard-scattering contributions to 5 — > P, y form 
factors. 

6.1 Low-energy representation of the B-type operators 

Below the intermediate scale jii ~ ^2EAqcd, the effective-theory description is in terms of 
soft, coUinear, and soft-coUinear messenger modes of SCETn [5, 14]. The Lagrangian in the 
soft sector is the same as in (5), while the coUinear Lagrangian becomes 

jOc^Cct(^^-Dc + ipc±Tj-f^ipc±]Cc- (100) 
2 \ m ■ JJc J 

Soft-coUinear messenger modes may be decoupled via field redefinitions in the leading-order 
soft and coUinear Lagrangians and in the four-quark operators representing the hard-scattering 
terms [15], and we do not display them here. In terms of the decoupled fields, it is convenient 
to introduce the gauge- invariant combinations [5, 28] 

Xc = [w^l ~ A , Qs = [s^Qs] ~ X'/' , :k = [S^h] ~ ^ 

yi^f - {iDii W,) ~ (A^ 0, A) , A^^ {iD^ S) ^ (0, A, A) , (101) 
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where Wc and S are coUinear and soft Wilson lines in the n and n directions, respectively. 
To account for non-localities generated by integrating out hard-collinear modes, coUinear and 
soft fields are smeared along light-like directions, i.e. (f)c{sn) and 0s (tn) [5]. 

Form factors arising in semileptonic and radiative B decays derive from current operators 
in the effective weak Hamiltonian mediating the decay of a heavy b quark into a left-handed 
light quark. The relevant operators are 

= g 7^(1 - 75) , Tr = i-t) qa^'il + 75) & = ? 7 (1 + 75) 6 . (102) 
We find it instructive when discussing radiative corrections to consider also the scalar current, 

SL = qil + lr,)b.^ (103) 

which can be related to the vector current by QCD equations of motion. In the NDR scheme 
with anti-commuting 75, the preceding analysis of operator mixing and rcnormalization is 
unchanged by the insertion of (1 + 75) next to the light quark. From now on, the replacement 
^ — > g (1 -I- 75) in the SCETj currents jf will be understood implicitly. The B-type current 
operators of SCETi then match onto four-quark operators in SCETn of the form [5, 6, 7, 8] 

Oi{s, t) = X,{sn) (1 + 75) Fc,, I Xc(0) Q,{tn) (1 ± 75) ^ F,,, %{0) . (104) 

The chirality of the field Qs{tn) is the same as that of Xc(0), which in turn is determined 
by the Dirac structure Fg^j. We restrict attention to the case where the matrix elements of 
these operators are evaluated with an initial-state pseudoscalar 5-meson. We then construct 
a basis of four-quark operators in which the soft component fields arise in the combination 
Qg (1±75) (^/2) D-C. (In a more general case, operators containing 7^^ next to Q,, and structures 
with more than one perpendicular Lorentz index, would also have to be considered.) For the 
scalar current Sl there is only one possible structure, 

05 = Xe(l + 75)|XeQ.(l + 75)|^H;, (105) 
while for the vector current we find the three operators 



o^, = ^x,(i + 75)^x, Q3(i + 75)|:k, 

O^y^ = t;^Xe(l+75)|Xe Q,(1 + 75)^:K, 

0(^3 = Xe (1 + 75) 7^^ I Xe Q. (1 - 75) I . (106) 

For the tensor current we only keep terms that are non-zero when contracted with q^,, where 
q with q±^i, = is the momentum transfer. This yields 

0^1= ^Xe(l + 75)5XeQ3(l+75)^3<, 
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0^2- ^5;e(l+75)7l'^X, Q,(1-75)^J{, 

n • V I I 

0^3 = v^^ % (1 + 75) I X, Q, (1 - 75) I J^ . (107) 

Color octet-octet operators have vanishing matrix elements between physical meson states, and 
also do not mix into the color singlet-singlet sector [15]. We may therefore restrict attention 
throughout to color singlet-singlet operators. 

In constructing the most general basis, we notice that adjacent 7_l matrices can always be 
avoided. Using 

{7^^,7l} = 2^r, (108) 

any such structure can be reduced to combinations of g^^ and totally anti-symmetric combi- 
nations. 

In four dimensions ^y^i - '^" vanishes for n > 3, while for n — 2 (we use e°^^^ = ~1) 

7M' = (f - f ) 75 ^ ^< (f - f ) 75 . (110) 

In dimensional regularization, so-called "evanescent" operators containing structures 
for n > 3 appear once radiative corrections are taken into account. A regularization scheme 
including the effects of these operators must be employed for next-to-leading order calculations. 
This is the two-dimensional analogue, in the space of transverse directions, of the procedure 
employed in four dimensions [29, 30, 31]. 

6.2 Matching calculations 

We introduce the momentum-space coefficient functions 

Di{E, cu, u,n) = J ds dt e-"' e'""'^-^ Di{s, t, n) , (111) 

where Di{s,t,n) is the coefficient multiplying the operator Oj(s, t) in (104), and the depen- 
dence of Dj on ri • P and n • v is only through the product (n ■ P){n ■ v) ~ 2E. The final 
step of matching SCETi onto SCETn is described by so-called jet functions Jij defined via 
the relation 

1 ( 2Eijj \ 

Di{E,uj,u,i^i) = J dyJij\^u,y, In— Cf{E,y,i^i) , (112) 

where are the coefficient functions of the SCETj currents J^, and iii ~ -y^EVVgcD denotes 
the intermediate matching scale. Since all interactions between coUinear fields and the b 
quark have been integrated out already in SCETj, the jet functions are independent of the 
heavy-quark velocity v. Dimensional analysis and rescaling invariance then imply that they 
depend on uj only through the dimensionless ratio 2Eu)/ iif. Furthermore, since only appears 
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logarithmically, the same must be true for this ratio of scales, to any finite order in perturbation 
theory [20]. This explains why above we have written the third argument of the jet functions 

as L = \n{2ELu/iJ,f). 

For the calculation of the jet functions, it is convenient to work in the primed basis of 
current eigenvectors jf defined in (42) and (44). From the structure of these currents, as well 
as the structure of the resulting four-quark operators in SCETn displayed in (105)-(107), we 
notice that 
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so that the SCETn operators Ovi,2, Oti and their SCETi counterparts Jv32, ^Ti ^^^^ related 
to the scalar operators 0$ and Jg simply by overall factors involving and v^. It follows 
that, e.g., the jet function arising in the matching of Jy^ onto Oyi is precisely the same as 
that arising in the matching of Jg onto Os- Similarly, we find that 
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so that the matching of J^j (Jtq) onto Ot2 (Ots) is precisely the same as that of Jy'^ onto 
Ov3- The operators Jyi and 7-^2 3 match onto SCETn four-quark operators with vanishing 
projection onto the i?-meson, while J-n 5 contain two perpendicular Lorentz indices. These 
operators are therefore not relevant to the B ^ P,V form factors, and we have thus not 
hsted the corresponding SCETn operators in (106) and (107). Collecting the elements Jij 
into matrices J', these observations imply that the jet functions can be expressed in terms of 
only two functions J'n and J'±. The result is 
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where the primes remind us that these results refer to the basis of the SCETi currents J^' . 
At tree level, we find 



J\\{U, V, L, /Xi)tree = J±{u, V, L, /ii)tree 



ATrCpasdii) 1 



N 



2Eu 



6{u — v) 



(116) 



In order to study the effect of next-to-leading order matching corrections, we now calculate 
the one-loop matching corrections to the jet functions J'\\ and J'±. The relevant diagrams are 
depicted in Figure 6. Contributions involving soft gluons cancel in the matching once the 
corresponding SCETn diagrams are evaluated, so that only hard-coUinear SCETi diagrams 
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Figure 6: SCETj graphs contributing to the matching of the currents J^' (crossed circle) onto 
the SCETii four-quark operators Oj. Full lines denote soft fields, dashed lines hard-coUinear 
fields. Diagrams with soft gluons or scaleless loops are not shown. 



need be considered. The 1/e poles are canceled by the renormahzation constants of the SCETi 
and SCETn operators. Details of this calculation are presented in [32]. The results take the 
form 



Jll{u,v,L,fii) 
J±{u,v,L,ij,i) 



N 2Eu 
iTrCpoisiiii) 1 



N 



2Eu 



-S{u-v) + ^^j^^{u,v,L) + 0{al] 
-S{u-v) + ^^j^{u,v,L) + Oial) 



where L — \a.{2Euj / ^f) , and the functions j\\ and j±_ are given by 

i\\{u,v,L) = ji{u,v,L) +j2{u,v,L) +j3{u,v) , 
j±{u, V, L) = ji{u, V, L) - j3{u, v) , 



;il7) 



(118) 
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with 
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The plus distribution is defined as in (48). At one-loop order there are also matching coef- 
ficients onto evanescent operators. We adopt a renormalization scheme in which the matrix 
elements of these operators between physical states vanish [29], and so do not list their match- 
ing coefficients here. However, the choice of evanescent operators is not unique, and the jet 
functions as well as the matrix elements of the physical SCETn four-quark operators depend 
on this choice. We have chosen the evanescent operators such that the matrix elements of the 
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physical operators are given by the MS LCD As of the B meson and hght meson [32].^ 

While the numerical impact of the one-loop matching corrections to the jet functions will 
be studied later, it is interesting at this point to use the explicit expressions in (118) to learn 
something about the "natural scale" to be used in the perturbative expansions for the jet 
functions. In the BLM scale-setting prescription, the scale in the couphng constant of the 
leading-order terms is chosen so as to absorb all terms proportional to Po — ^ Ca — ^TpUf 



in the next-to-leading order corrections 
and uj values) we should replace as{i-ii) 



In the present case, this means that (for fixed u 
Q!s(a*blm), where 



/^Ilm = 63. 2uEuj 



(120) 



for both jet functions, which agrees with previous results reported in [34]. When applied to a 
specific process, such as the hard-scattering contributions to heavy-to-light form factors to be 
considered in Section 7, the BLM scale setting should be done after the convolution integrals 
over the LCDAs have been performed, which is equivalent to using mean values u — > e^'""^ 
and u g(ina;)^ weighted by the corresponding tree- level amplitudes, to define an effective 
scale for each decay process. The resulting BLM scales, while parametrically of order /Xj, 
are numerically rather low, typically /xblm ~ 0.7 GeV or less. The smallness of these scales 
appears to suggest a poor convergence of the perturbative expansions of the jet functions. 
However, it turns out that the BLM prescription cannot be trusted in the present case. In the 
limit where we neglect the mild y-dependence of the coefficients {E,y,fii), it follows from 
(112) and (117) that the convolution integrals relevant to the hard-scattering contributions 
governed by the longitudinal jet function J'\\ yield 
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where (...) denotes an average over the 5-meson LCDA (f)B{uj,fii) with measure du/uj, and 
for simplicity we have assumed the asymptotic form (f)M{u, /ij) = Quu for the LCDA of the light 
meson. The corresponding result for the jet function J7j_ is obtained by replacing jy j^, in 

2 

which case the coefficient of the single logarithm changes to (—6-1- ^), and the constant term 
changes to 0.73. If scale /li is chosen such that the logarithms appearing in this expression 
are not too large, then (at least through one-loop order) the perturbative corrections are of 
moderate size. For comparison, keeping only the BLM terms of order Pooi^, we would obtain 



«£(/fi) 

An 



25 /, 2Euj\ 
— ( In ) + 26.39 
3 \ /^i / 



(122) 



instead of the expression in the second line of (121). We conclude that the BLM prescription 
overestimates the size of the 0{a'^) corrections by large amounts. The perturbative expansions 
of the jet functions at a scale //, = \/2EAh^ 1.5 GeV appear to be reasonably well behaved. 

^We are indebted to M. Bcnckc and D. s. Yang for pointing out that this was not true with the choice of 
evanescent operators adopted in the preprint version of this paper. 
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6.3 Sudakov resummation in SCETn 



Having performed the matching onto SCETn, we now wish to evolve the momentum-space 
coefficient functions A in (111) from the intermediate scale jii ~ ■yj2Eh.qcT) down to a 
hadronic scale yU, which is independent of the large energy E. This can be achieved by solving 
the integro-differential RG equation [15] 



— Di{E,u!' ,u , = / du du^r{i^,^^' ,u,u' , fi) Di{E,uj,u, 



(123) 



The anomalous dimension 7r depends on the spin-parity quantum numbers of the light final- 
state meson M . As before, we must distinguish two cases: T = || for M — P,V\\, and F = ± 
for M = Vj_. The anomalous dimension may be decomposed as 



where 



7r((^, oj' I u, u', /i) — S{u! — cu') Kr{u, u') -\- S{u — u') H{u!, cu', /i) , 



vv Ki'{u, v) — 



TT 



_ 9(v — u) _ 9(u — v) 
uv h vu 



V — u 



u — V 



H — uuS{u — v) 
2 



(124) 



- cr 



uve{v - u) +vu9{u- v) \ + 0{a^) , 



with C|| = 1, cj_ = is the Brodsky-Lepage kernel [24], and 

H{u, uj\ jj) = rcusp(Q;s) In — j{as) S{lu - lu') + lu T{uj, uj', as 

L CU -1 



TT 



u{u — uj') uj'{uj' — uj) 



(125) 

(126) 
+ 0{al) 



is the corresponding kernel governing the evolution of the 5-meson LCDA [20, 36]. 

As a consequence of the conformal symmetry of QCD at the classical level, the kernel 
vv Ky{u,v) is symmetric in u and v, and Ky is diagonalized by Gegenbauer polynomials (cf. 
(96) in Section 5.3). We write 



/ duKY{u,v)ipn{u)^ KY,nVn{v) , 

Jo 



where at one-loop order 



TT 



mn^l - 1) + ^ - + 2) 



(127) 



(128) 



We choose the normalization of ipn such that they are orthonormal in the measure uu du on 
the unit interval. 



2n + 3 



(n+l)(n + 2) 



CrH2.-l) = y^^^^^p(^'^)(2.-l). (129) 
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After expanding the coefficients Di{E,u!,u, /i) on the basis of eigenfunctions (pn{u), the RG 
equation for each Gegenbauer moment is solved as in [20]. The result reads 



Di{E,uj,u,ii) = 



2Euj \ii 



(130) 



n=0 



X 



/ dyCf'{E,y,iii) / dx xx (fin{x) Jlj{x,yyr„ l^'i 
Jo Jo 



2EluV 



7?=0 



where S{ii, /li) and a{fj,i, /i) are given in (64) and (65), the jet functions J^lj in the primed basis 
are obtained from (112) by replacing 2£'a;//x? — > V,,, and 



Vr,n{f^i,H,rj) 



J^{r] - a{iii, ii{a)), a) + KT,n{cy) + 7(«) 



(131) 



with 



J^{r],as) = J dujuj'r{LLi,uj',as) (^^^ 

= ri-5°P(a,) [^(1 + rj)+ ^'(l -r)) + 2^ + 0(al) . 



(132) 



For a complete leading-order solution, we require the two-loop expression for Fcusp appear- 
ing in S{fj,,fii), one-loop expressions for the remaining anomalous dimensions in a{iJ,i, /i) and 
Vr,n{l^i, l^,v)j tree- level matching conditions for the jet functions. At this order, the 
relevant expansions are 



4A^ 



^r,n(A*i,A*, 0) = In 





To 



4:71 

lnr2 , 



:i-r2 + lnr.) + AlnV,-(^-|) (l-l + ln 
2/?o \To poj \r2 



r2 



2Po 



y-7i.^lnr,-^(4°i + 7o)lnr,. 



(133) 



where r2 = as{iJ>) / OisifJ'i) ^ 1; and as usual we have expanded 



(0) "s , , 

4,L^ + --- , 7 = 70^;; + .... 



47r 



(134) 



7 Application to heavy-to- light form factors 



In this section we apply our results to obtain resummed expressions for the hard-scattering 

contributions to heavy-to-light form factors at large recoil energy. We begin by recalling the 
definitions of the ten form factors describing B decays into pseudoscalar and vector mesons. 



35 



Equating these definitions to the corresponding effective-theory expressions shows that at 
leading order in 1/E only eight of the form factors are independent. As shown in (1), in this 
approximation the form factors are given by the sum of a soft-overlap contribution, expressed in 
terms of universal non-perturbative matrix elements (m{E) (where M = P, V\\, V± depending 
on the spin-parity quantum numbers of the light final-state meson), and a hard-scattering 
contribution, given by a convolution integral over meson LCDAs. We focus on the second term 
and present general, resummed expressions for the corresponding hard-scattering kernels. The 
results are particularly simple when hard matching corrections arc ignored. In this case, for 
a given final-state meson M, the hard-scattering contributions to all B ^ M form factors 
can be parameterized in terms of a universal quantity Hm- The Sudakov suppression factors 
are mild in all cases. In order to investigate the relative importance of running and matching 
corrections, we apply our previous one-loop matching results for the scalar current to obtain 
the one- loop corrected hard-scattering contribution to the vector-current form factor Fq, and 
to the form factor at maximum recoil {q^ = 0). 



7.1 Form- factor definitions and factorization formulae 



We begin by recalling the definitions of the form factors parameterizing B decays into pscu- 
doscalar (P) and vector {V) mesons, following the conventions of [13] (we use the sign con- 
vention e°^^^ = — 1): 



{P{p')\q7'b\B{p)) = Mq') I P'+P'' - 



m%-mp 



^2, ml -ml 



2 qn+Fo{q') "2 ' 



{P{p')\qa^''q.b\B{p)) 
{V{j)\r))\qYh\B{p)) 



mB + mp 

2iV{q^) 
mB + my 



[qHp'+p'n-{ml-ml)q^], 



e^'^'W.p'pP., 



{V(p', v)\qri5 b\B{p)) = 2mvAo{q') q" + (mB + my) A^{q') (77*'^ -"^q^ 



A2{q' 



q^ 

7]* ■ q 



ms + my 

{V{p\r^)\qa^^q,h\B{p)) = -2T,{q^) e^^^'^r^lp'p, ^ 



q^ 



{V(p',r))\qa^^^,q,b\B(p)) = -iT,(q') [(ml - m'y) rj*^ - rj* ■ q (p'^ + p'^)] 



iT^{q )ri* -q 



q' 



mo — m 



{]f + p'>') 



V 



(135) 



Here q = p — p' is the momentum transfer, and 77 denotes the polarization vector of the vector 
meson, which satisfies rj • p' — and 77-77* = —1. The polarization vector for longitudinally 
polarized vector mesons is given by 



P 



my n ■ p' 



7?T.y n ■ V 



(136) 
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For the evaluation of the hadronic matrix elements entering the hard-scattering contribu- 
tions to the form factors, we need the definitions of the leading-twist LCDAs for the light 
mesons in terms of SCETn fields. They are [9] 

(P(p')|Xc(sn)r|Xe(0)|0) = ^n-p' tr (^^^T^^ e^^^'^^P , 

{V^ip', ri)\X,{sn) r I X,(0)|0) = '-^^ n ■ p' tr (^f^ T ^'dx^e-^^'^' cf>vAu, /.) , (137) 

where the 0m are normalized according to J^du(l)M{u, fi) = 1. Note that the transverse 
vector- meson decay constants are scale dependent. Below, we will sometimes write /A/f(/x) as a 
generic notation for the light-meson decay constants, keeping in mind that a scale dependence 
is present only in the case of fvj^- The LCDA for the S-meson is defined as [35, 36] 

{Omn)^rnom = -'^^/^^M tr ri±^75) j^du;e-^-'-^M^.lA , (138) 

where denotes the S-meson state defined in HQET, and F is related to the asymptotic 
value of the 5-meson decay constant in the heavy-quark limit (see below). 

For completeness, we list also the soft-overlap contributions to form factors, parameterized 
in terms of universal matrix elements (m [12, 13]. At leading order in 1/E, we define^ 

{P{p')\X,,T h\B,) = 2ECp{E) tr (^^rii^) , 
m\ip')\X,,T h\B,) = -2ECv^^{E) tr (^^ri±^75) , 

{V^ip')\X,,T h\B,) = 2ECvAE) tr (^C^r^^s) • (139) 

The definitions are given in terms of SCETi currents. As shown in [8], the functions Cm{E) 
can be decomposed further into matrix elements of SCETn operators. These operators all 
satisfy the same symmetry relations, and the linear combinations contributing to the form 
factor behave under renormalization in precisely the same way as the SCETj operators used 
in the definitions (139). 

To simplify the factorization formulae, it proves convenient to write the coefficient functions 
in the form 

Di{E, u, u,i^) = - ^^I^ Ti{E, u, u, li) , (140) 

^ mBjB 

^Our convention for (v^^ is the same as in [12], and is such that all Cm functions have the same power 
counting. This differs from the corresponding function ^|| in [13], given by Cvy = (E/mv) ^\\- 
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which serves as a definition of the hard-scattering kernels Tj. Here Kp relates the HQET 
parameter F in (138) to the physical i?-meson decay constant, fB^jT^^B — Kp^ij) F^ji), up to 
corrections of order h.qcT>lfT^b- At next-to- leading order 



Kf{ij^) = 1 + 



An 



31n^-2 



;i4i) 



With all of the definitions in place, it is now a straightforward matter to equate QCD 
matrix elements to the corresponding SCET expressions. For the scalar and pseudoscalar 
currents, for example, we obtain 

— {P{p')\qb\B{p)) = —Ci{E,i,)Cp{E,i,) 
ruB mB 



+ 



IE dcj 

/ — (I>b{(^,I^) / dufp(t)p{u,jji)Ts{E,uj,u,ij) 

rriB Jo Jo 

— [CiCp + 'l>B^fpcl>p^Ts] , 
niB 



1 - 9 F" r 



Vn 9Vn 



Tc 



(142) 



where in the second line we have introduced a symbolic notation for the convolution integrals. 
In general, the scale /i in the soft-overlap terms could be taken different from that in the 
hard-scattering terms. Similarly, from the matrix elements of vector and tensor currents we 
find the relations 
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f 1 - — ) Tt3 




V mBj 



mB 
2E 



To 
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Tl 
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2^ 



ruB J 



T2 



a 



T3 



+ ( 1 ) 

TUB J 



mBT2-2ET3 
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/^A ^ /^A I -'- /^A 



r4 



/l 



Tl 



(143) 



We have dropped terms of relative order {vriM/EY-, where mjvf is the hght meson mass. Recall 
that we have defined E = v ■ P^, which differs from the true energy by an amount of order 
m\j/E. Also, at leading order in Aqcd/'"^6 there is no difference between and mf,. The 
results for F+, Fq and in (143) agree at tree level with the corresponding expressions in 
[19].^ Beyond tree level, and to any order in perturbation theory, we find that there is still only 
a single jet function for decays into a given light final-state meson, e.g. J\\ for a pseudoscalar 
meson. This implies that the two functions Ja and Jb introduced in [19] are not independent. 
Our result has important implications for the universality of hard-scattering corrections, which 
we discuss in more detail in the following subsection. The general results (143) establish that 
the form-factor relations 

A(g^) = r^^^ V{q-) , Uq-) = ^ T,{q-) , (144) 
(mB + mv) TUB 

are rigorous predictions of QCD at leading order in 1/E, and to all orders in Ug. These 
relations were conjectured in [39], and they were shown to hold at one-loop order in [13]. 

The scalar and pseudoscalar currents may be related to the vector and axial- vector currents 
using the equation of motion for the quark fields, as in (31). When applied to the matrix 
elements appearing in (135) and (142), it follows that (setting as previously the light quark 
masses to zero) 

^fo(/^QCD) 



ml 



(P(p')|[g-&](/^QCD)|S(p)) = 1 



nil 



nil 



Fo{q')^Fo(q') + ..., 



m. 



(M|(p')|[?75&](/^qcd)|S(p)) = 



2m 



2E 



^(?2) + ..., (145) 



niB 

where the dots represent terms of order m\j/m\. Comparison with (143) shows that to all 
orders in perturbation theory 

2 I I ^A 

nib) 



'V2 



+ a 
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V3 



Tvi + 



1 - —) Tv2 = 

nib J 



nib 

rfib{iiqcT>) 
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s > 



(146) 



These are the SCETn realizations of the SCETi coefficient relations found in (31). 

[38]. 



'See also [37]. In equation (22) of this paper 



-^4*^ should be replaced by +B^^^ 



39 



7.2 Sudakov resummation for the hard-scattering contributions 

We now apply the results of Sections 5 and 6 to obtain completely general, resummed expres- 
sions for the hard-scattering kernels Ti appearing in the form- factors (143). Starting from the 
defining relation (140) for Tj, and combining it with the result (130) for the coefficients A, 
we obtain 



Ti{E,uj,u,^) = ————— [ — - e 

8E KpiiJih) Kpt-hJ VA* 



oo 

X 

n=0 



X 



oo „i „i 

^ipn{u) dy dvUT{y,v,iJ,h,IJ'i)Cf{E,v,iJ,h) 
„=o Jo Jo 

\ J 



(147) 

»7=0 



We have used (63) to express the coefficients C^' {E.y, /ii) at the intermediate scale, which 
enter in (130), in terms of their matching conditions at the high scale Hh- Also, we have 
rewritten Kf{h) = e"^^^'^'"'^) Kpinh), where 

VFi^in, //) = - r^'^ 7f(«) = ^ Innr, + ... (148) 
is the solution to the evolution equation 

7 °0 1 1 

Kf = ^fKf, 7F = J]7F,n(^j , (149) 



dluji 



n=0 



with ^Ffi — — 3Ci?. As previously, ri = as{i^i) / (y.s{nh) and r2 = as{n) / as{ni) . In (147) all 
large logarithms are resummed into the various evolution functions. The matching coefficients 
Kpifih), Cf {E,v,iih), and J'-j{x,y,'Vr^, fii) can be calculated reliably in fixed-order pertur- 
bation theory. The expressions for the kernels Tj are formally (to the order we are working) 
independent of the two matching scales and /Xj. A dependence on the low scale /j, remains, 
which will cancel against the scale dependence of the hadronic matrix elements of the SCETn 
operators. 

The general result (147) simplifies considerably when hard matching corrections at the 
high scale fj,h ~ nib are ignored. In this case Kp{fih) = 1, and using the fact that the tree-level 
matching conditions for {E, v, Hh) are independent of v, it follows that 

»i 

dvUr{y,v,iih:l^i)Cf'{E,v,iih) = Ur{y, i^h: l^i) Cf {E, ^h) , (150) 



where Ur{y, Hh, l-^i) are the functions shown in Figure 4. Using the tree-level matching con- 
ditions compiled in (26)-(28) along with the relations (43) and (45), it is straightforward to 
evaluate the products J'lj Cf {E, fih) for the various kernels Tj, with the matrices Jl^ as given 
in (115). We find 

Ts- -J\\, 
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Tvi ■ ( 1 ] '^V2 ■ - 2 Jji , Tv3 ■ , 

V mbj 

AE 

Tti : 2J|| , Tt2: —J±, Tts : . (151) 

rrib 

It follows that the convolution integrals for the hard-scattering contributions to the form 
factors in (143) may be written in terms of three universal functions 



Hm^-^( — ] e^('^-'^^)+^('^''^')+^-('^-'^) V / du fMifi) <Pm{u, fi) 

X / — 0B(a;,/x) - / dyUr{y,Hh:f^i) 

Jo ^ V/^/ ^0 



/ dx XX ipr,{x) Mx, y, Vr,, l^i) e^^-^^^^''"'''^ 
Jo \ f^i / 



(152) 

■n=o 



where T = \\ for M = P or V\\, and F = ± for M = V±. We have extracted a factor —{itlb/IEY 
from the expression in (147) so as to remove any non-logarithmic E dependence from the 
quantities Hm and ensure that they are positive. (Recall that the tree-level jet functions are 
proportional to —1/E.) Note that in (152) the ji dependence of the kernels Tj cancels against 
that of the hadronic quantities 0s, 0m, and /m- We repeat that this expression neglects hard 
matching corrections at the scale ~ but it allows for arbitrary corrections to the jet 
functions. If only the tree-level matching conditions for the jet functions given in (116) are 
retained, the result simplifies further to 

/I I y Li'i r ^ 1 I f,~' 1 I LJ I / . I ■ J \ 



N mB\fihJ f^Jo 

^ ^v,M+Vr,U^,,^,o) dyy^r^{y)Ur{y,|Ih:l^^). (153) 

Jo ^ V/^/ Jo 

Prom (143) and (151), it is easy to read off the hard-scattering contributions AFj to the 
various form factors. For the three S — > P form factors, we obtain 



The corresponding results for the seven B —> V form factors read 



2m 

AA2 = 



TUB 
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and 

^^^^ 2^l^^- + ^^^^ij- ^^^^^ 
7.3 Numerical results 

We are finally in a position to study the phenomenological implications of our results. For the 
leading-twist LCDAs of the light mesons, we take for simplicity the asymptotic forms 

(f)p{u) = {u) = ("w) = 6^(1 - u) . (157) 

Then only the n — term contributes to the sum (152), where ipo{u) — y/6. To proceed 
further we require the S-meson LCDA, which is poorly known at present. As an illustrative 
model, we adopt the form [35] 



M^, /^o) = ^ e-"/^^ , XB = l{mB- nib) ~ 0.32 GeV , (158) 



where Hq is a low hadronic scale, at which the model assumes the stated functional form. The 
relevant convolution integral resulting from (153) at /i — /iq is then 



—a 



^M-,^o) (-) =f r(l-a)(^) . (159) 



The quantity 



Combining all pieces, we find at leading order in RG-improved perturbation theory 

N 1713X3 \l^h) V/^o, 

,v^{/..,^.o)+^o(w,^o) j\y 2y Ur{y, i^h, l^i) . (160) 
Jo 

V„M = l„r(l - I| m..) - (-,.5j + J£) inr.. (161) 

where now r2 = as{fJ^o)/as{fJ'i), results from the product of the F function in (159) times e^'"-" 
from (133), after the term involving has been used to cancel the scale dependence of 
/m(a*), replacing it with fnilJ'i)- We expand the RG factors a{iih, /-li) and S{iih,IJ'i) according 
to (66), a(/ii,/io) and S{iio,iii) according to (133), and VF{f^h, f^) as shown in (148). 

It is convenient to present the results in terms of the tree-level expressions, which we define 
at a fixed reference scale = y/rribAh with A/^ = 0.5 GeV, 
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Figure 7: Leading-order resummed results for the hard-scattering contributions to heavy-to- 
hght form factors evaluated at 2E = m^, in units of the tree-level expressions for H]^^^ in (162). 
The bands reflect the sensitivity to the choice of the two matching scales Hh (dark band) and 
Hi (hght band). They are the result of varying 0.5ml < A*fe < and O.bnibAh < l^i < 2mbA/i, 
respectively. 



As an example we may consider the pion (with = 131 MeV) and the p meson (with = 
198 McV and fp^ (Vm^A/i) = 152 McV [40]) as representative pseudoscalar and vector mesons, 
respectively. Then with fe = 180 MeV and = 0.32 GeV, we obtain ^ 0.021, H^^^^"" ^ 
0.032, and Hy^^ 0.025. In Figure 7 we plot the dependence of the quantities Hm on the 
model parameter /iq and investigate their sensitivity to the matching scales Hh and //j. We 
observe that the results are rather stable for a wide range of /xq values. If 0s can be modeled 
accurately by (158) for some value of Ho within this range, then the results are relatively 
insensitive to the precise value of hq. At /xq = 1 GeV, the sensitivity to the matching scales is 
approximately ±20% for the considered range of /^j, and ±5% for the considered range of 

In order to investigate the size of one-loop matching corrections, we return to the example 
of the scalar current. Inserting the one-loop jet functions calculated in Section 6.2 into the 
general relation (147), it is straightforward to project out the relevant moment of the light- 
meson LCDA. The second relation in (146) then yields the hard-scattering contribution to the 
form factor Fq including leading-order resummation effects and next-to-leading order matching 
corrections. Figure 8 shows the energy dependence of AFq for two different values for the model 
parameter hq. In Figure 9 we restrict attention to maximum recoil energy, E = mB/2. In 
this case -^+(0) = Fo{0), so that our analysis applies to both vector form factors. Wc again 
observe a mild dependence on the model parameter /iq, and a slightly reduced (with respect 
to the leading-order results in Figure 7) sensitivity to the matching scales i^h and /Xj. 

As a final note, we emphasize that the numerical results presented above depend on the 
model used for the S-meson LCDA, primarily through its first inverse moment, A^. However, 
the statement that leading-order resummation effects are encoded by universal functions Hm, 
which are the same for all form factors with the same light meson in the final state, is inde- 
pendent of any model assumptions. The one-loop matching analysis for the sample case of 
the form factor Fq provides some insight into the size of next-to-leading order corrections and 
the convergence of the perturbative expansion of the jet functions. The one-loop matching 
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Figure 8: Hard-scattering contribution to the form factor Fq including one-loop matching 
correction, in units of the tree- level expression AFq'''^'' = ^ Hp'^'^. The solid and dashed 
curves correspond to a low-energy scale /iq such that a;<i(/xo) = 0.5 and 1.0, respectively. The 
black hues are obtained with matching scales Hh — 2E and //^ = y/2EA^ with = 0.5 GeV, 
the gray lines with /i^ — irib and //j = \^mbAh. 
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Figure 9: Hard-scattering contribution to the form factor F+(0) = Fq{0), including one-loop 

matching corrections, as a function of the model parameter ^q. The bands reflect the sensitivity 
to the choice of the two matching scales fih (dark band) and fii (light band). They are the 
result of varying 0.5m^ < ii\ < 2ml ^^"^ O.SmftA/j < < 2mfeA/i, respectively. 
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corrections at the hard scale modify the tree-level expression for Fo(0) by ~ 10%. With 

/io = 1 GeV, the corrections arising from matching at the intermediate scale fii amount to a 
change of ~ 30%. Our calculation also confirms the convergence of the hard-scattering con- 
volution integral in a non-trivial required by the factorization theorem (1), and as 
proved on general grounds in [7, 8]. 



8 Asymptotic behavior of Sudakov logarithms for 
heavy-to-hght form factors 

An interesting question to ask is whether one of the soft-overlap or hard-scattering contri- 
butions to the form factors is suppressed relative to the other in the formal limit E ^ oo 
with fixed ratio E/rrih. As already noted in [8, 15], this issue cannot be entirely addressed 
in perturbation theory, since the matrix elements for the soft-overlap contribution depend on 
in a non-perturbative way. Still, it is interesting to investigate the relative suppression of 
the Wilson coefficients multiplying the hadronic matrix elements at some low renormalization 
scale /io, assuming that the matrix elements at a low scale have only a mild energy dependence. 
In that way, one takes into account all short-distance logarithms arising from RG evolution 
between the scales /ih ^ and /^o- 

In order to isolate the dominant Sudakov factor, we concentrate on only those terms in 
the RG evolution equations involving Fcusp- They are sufficient to resum the leading Sudakov 
double logarithms to all orders in perturbation theory. In SCETi, the evolution of the A-type 
and B-type operators is then identical. 



rcusp(«s) In 



2E 
1^ . 



Fcusp In — 



Cf{E,ii), 



Cf{E,u,i,), 



(163) 



so that the leading Sudakov suppression factor from SCETj running is the same for the two 
types of currents. Our new result (54) for the cusp contribution to the anomalous dimensions 
of the B-type currents is a crucial ingredient to this conclusion. In SCETn, the leading-order 
heavy-light currents contributing to the soft-overlap terms obey the same evolution equation as 
in (163) [8]; however, the hard-scattering contributions now derive from four-quark operators, 
whose coefficient functions obey 



— ^A(-E, a;, !*,//) ~ [rcusp(Q;s) In - Di{E,u;,u, /j) 
a in /i L 00- 



rcusp(as) In 



2E 
1^ 



rcusp(Q;s) In 



2Euj' 
1^. 



Di{E,uj,u,iJ,). (164) 



If not for the term involving the intermediate scale 2Eu; ~ E'Aqcd ~ A*f , the running would 

again be the same for both contributions. For a low value /i ~ /iq ~ Aqcdi the effect of this 
term is to reduce the Sudakov suppression of the hard-scattering contributions relative to that 
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Figure 10: Asymptotic Sudakov suppression of the Wilson coefficient of the soft-overlap con- 
tribution to a heavy-to-light form factors relative to the coefficient of the hard-scattering 
contribution. The three curves correspond to different choices of the low scale /^o, such that 
Q!s(a*o) — 0.5 (solid), 0.75 (dashed), and 1.0 (dotted). The intermediate matching scale is 
taken to be jii = ^/2Ekh with = 0.5 GeV. 



of the soft-overlap terms. It follows that, at a low scale <^ ^i, the soft-overlap contribution 
to a heavy-to-light form factor is suppressed relative to the hard-scattering contribution by a 
factor 

""=('^0) da r da' 



exp 



-2 



_,25(Mi,/io) 



(165) 



where the leading-order expression for S{iii, hq) is obtained from (66) by replacing fj^h fJ^i 
and II ^ /Iq. For realistic values of parameters this suppression is very mild. This is illustrated 
in Figure 10, where we show the energy dependence of e'^^^^^'^^ for /ij = y/2EAh and different 
choices of /Iq. 



9 Discussion and conclusions 

Our analysis completes the leading-order description of large-recoil heavy-to-light form fac- 
tors in SCET, including Sudakov resummation effects. The main focus was on the struc- 
ture of SCETi currents, which divide naturally into two classes (denoted A- and 5- type), 
contributing respectively to the soft-overlap (spin-symmetry preserving) and hard-scattering 
(spin-symmetry breaking) terms. 

Two interesting complications arise in the RG analysis in SCETj. The ffist, common 
to both types of operators, is the explicit scale dependence of the anomalous dimensions. 
The non-trivial result that the scale dependence appears to all orders in as only as a single 
logarithm times the universal cusp anomalous dimension follows from the cusp structure of the 
Wilson loops appearing in the effective-theory operators. The resulting Sudakov suppression 
is the (formally) dominant effect of RG running. The second complication, specific to the B- 
type operators, is the mixing of operators having different hard-coUinear momentum fractions. 
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The structure of this mixing is constrained at leading order by a conformal symmetry in the 
hard-colhnear sector of the effective theory. The corresponding (multiphcatively renormahzed) 
operator eigenfunctions are expressed in terms of Jacobi polynomials. The heavy quark breaks 
the conformal symmetry, introducing an additional local term in the evolution equation. An 
approximate, but arbitrarily precise solution to the RG equation may be obtained by matrix 
methods with a truncated basis of conformal eigenfunctions. 

To complete the form-factor analysis, we have matched the intermediate effective theory 
onto SCETii, where the i?-type currents are represented in terms of non-local four-quark op- 
erators. The projection properties of effective-theory spinor fields severely restrict the possible 
Dirac structures appearing in the four-quark operators. As an immediate consequence of ex- 
pressing the original QCD currents in the SCETn representation, we found that for two cases 
the spin-symmetry relations holding between the soft-overlap contributions to different form 
factors are not broken at any order in perturbation theory by hard-scattering corrections. In 
the remaining cases, when matching corrections at the high scale /ih ~ are neglected, the 
resulting hard-scattering contributions can be related to a single factor Hm, the same for all 
B ^ M transitions involving the same light final-state meson M. This universahty imphes 
that certain form-factor relations, e.g. those between the B ^ P form factors F^, Fq, and F^, 
are exact up to corrections of order as{mb) and Aqcd/^6- We have presented results for the 
perturbative expansion of the quantities Hm through one-loop order. 

Because of the factorized form, the RG analysis of the hard-scattering terms in SCETn 
reduces to the renormahzation of the 5-meson and hght-meson LCDAs. Using previous results 
for the relevant evolution equations, we have presented a complete leading-order resummation 
of Sudakov logarithms. At this order, after resumming all single and double logarithms, 
the result is described by a universal RG factor, identical for all form factors describing 
the same light final-state meson. In order to investigate the size of matching corrections, 
we have performed the one-loop matching necessary to analyze the vector form factors F+ 
and Fq at maximum recoil. Contributions from matching at the high scale fih ~ 2£' ~ 
give ~ 10% corrections to the tree-level expressions, while matching contributions at the 
intermediate scale /li ~ ^2£^Aqcd yield modifications of ~ 30%. A complete next-to-leading 
order solution, which controls scale dependence through 0{as), will have to await the extension 
of our anomalous-dimension calculations to the two-loop order. Such an analysis would also 
require the three-loop coefficient of the cusp anomalous dimension rcusp(Q;s), which has been 
calculated very recently [41]. 

One of the dominant uncertainties in the phenomenological analysis of heavy-to-light form 
factors results from our ignorance about the functional form of the S-meson LCDA 0b (o;). 
The scale of the hard-scattering contribution to form factors is set by the quantities Hm ~ 
(47r/3)ajB/M/mijAs ^ 0.02-0.03 (for a, ^ 0.3 and ^ 0.3GcV). Since (puito) and, in 
particular, its first inverse moment are ubiquitous in the effective-theory description of 
exclusive B decays, there is hope that experiment can provide significant constraints, to be 
checked against the present knowledge based on QCD sum rules [35, 42, 43]. 

The effective theory does not predict the relative size of the soft-overlap and hard-scattering 
contributions to the form factors, other than saying that they are of the same order in Aqcb/E 
(neglecting Sudakov logarithms), namely of order (Aqcd/-^)^^^- To explore this question 
further, we have investigated the asymptotic forms of the Wilson coefficient functions in the 
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formal limit £^ — > oo with fixed E/rrii,. We have found that the soft-overlap coefficients are 
suppressed relative to the hard-scattering coefficients in the asymptotic large-energy limit, 
but that this suppression is ineffective for realistic energies attainable in B decays. This 
conclusion is, however, tempered by the fact that the non-perturbative soft-overlap matrix 
elements have a long-distance sensitivity to the scale E, which cannot be controlled using 
short-distance methods. Therefore, the behavior of the Wilson coefficients alone does not 
necessarily give a reliable prediction for the asymptotic behavior of the form factors themselves. 
Phenomenologically, it appears that the soft-overlap contributions to B ^ M form factors are 
significantly larger than the hard-scattering terms. 

The results presented in this paper will be relevant to more complicated decay processes. 
QCD factorization formulae relate the decay amplitudes for rare exclusive processes such 
as S — > 7171 or B ^ K*'j to the sum of a S — > M form-factor term plus a hard-scattering 
contribution [9, 10, 11], both of which are described in the effective theory by operators already 
present in the form-factor analysis. For example, in [44] it is argued that the jet functions 
appearing in S ^ tttt are the same as those appearing in the S — > tt form factor. The 
universality of the jet functions discussed in the present work may have interesting implications 
when applied to these cases. Finally, using the SCET formalism it should also be possible to 
analyze heavy-to-light form factors beyond the leading order in the large-energy expansion. 

Note added: While this paper was in writing the work [45] appeared, in which the one-loop 
hard matching corrections to the subleading SCETi currents are computed. Our operator 
basis is more convenient for studying RG evolution than the one adopted in that paper, since 
with our choice the two-particle (A-type) and three-particle (i?-type) currents do not mix 
under renormalization. Our matching coefficient Cg in (30) agrees with the corresponding 
result in [45]. 
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